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SUMMARY 


The electronic controller described in this report, together with an 
active control landing gear, are found to provide significant reductions in 
forces sustained by an aircraft during takeoff, landing impact and rollout. 
The degree of force reduction increases with the sink rate. These results 
were obtained analytically and confirmed experimentally by actual drop tests 
of a landing gear under active control. 

The electronic controller continuously compares the kinetic energy of 
the aircraft with the work potential of the gear until the work potential 
exceeds the kinetic energy. The wing/ gear interface force present at this 
condition becomes the command force to a servo loop which maintains the 
wing/ gear interface force at this level by providing a signal to an electro- 
hydraulic servovalve to port flow into or out of the landing gear. 


INTRODUCTION 

Hydraulic Research (HR) was retained under NASA Contract NAS 1-14459 
to design, develop, fabricate and test an electronic controller for an electro- 
hydraulic active control landing gear, as described in Reference 1. The primary 
function of the active control landing gear is to minimize the force to which the 
aircraft is subjected as a result of landing impact and rollout, takeoffs and taxi 
operations. As shown in Ref. 1, the resultant decrease in applied loads could 



be very beneficial in reducing structural fatigue design problems and increas- 
ing ride comfort for passengers and crew. The work was divided into two 
major phases: Phase 1 for the analytical development and design of the 
controller; and. Phase 2 for the fabrication and testing. This report summa- 
rizes the effort of both phases and provides the rationale for the controller 
design selected. 

Phase 1 included the development and use of digital computer programs 
which simulated the dynamics of the landing gear/ controller system. A non- 
linear dynamic model and a simplified linear model for vertical drops were 
developed by HR and used extensively in the analysis. In addition, a computer 
program developed by NASA, which included aerodynamic simulation and 
landing gear dynamics, was supplied to HR in support of the program. The 
computer program was modified by HR to include the dynamics of the electro- 
hydraulic servovalve and this program was incorporated into the study. The 
landing gear parameters, as described in the NASA computer program, were 
used as the basis for all analyses described herein. During this phase, the 
controller design was also developed. This design included analog circuitry 
for the control laws, compensation, summations and the signal to the servo- 
valve driver while a digital microprocessor was used for the nonlinear 
computations and decision making. 

During Phase 2, the controller, as defined in Phase 1, was fabricated 
and tested by actual drop tests of the gear. 
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SYMBOLS 




A 0 area of orifice in shock strut orifice plate, 0.0000786 m a (0. 1219 in 2 ) 

A! shock strut hydraulic area (piston area), 0.00317 m a (4.909 in 2 ) 

A s shock strut pneumatic area (cylinder area), 0.00535 m» (8. 286 in 2 ) 

A 3 annular area in shock strut between piston and cylinder walls, 0. 00151 m a 

(2. 34 in®) 

discharge coefficient for active control servovalve orifice, 0.62 

discharge coefficient for shock strut orifice, 0. 90 

C Q orifice coefficient for shock strut orifice = Cj 0 A 0 V 2g c / p , 3.45 (10" e ) 

m 4 sec' 1 • ( 17.51 in 3 /sec/ V lb/in a ) 

CP linearized orifice coefficient for active control servovalve ='9Q SV /3P 1 

= 3. 16 (10‘ 1X ) m 6 • N' 1 • sec' 1 (0. 01334 in 3 / sec /lb /in 2 ) 

CP G linearized orifice coefficient for shock strut orifice = 9Q 0 /SP 1 
= C 0 / (2 P 1 -P a ), m B • N’ 1 • sec' 1 (in 3 /sec /lb/ in 2 ) 

CQ linearized orifice coefficient for active control servovalve = SQ sv /dX sv 
^sv V(P S + P r )/ 2, 8.61 m a / sec (13 340 • in 3 /sec/in) 

C sv orifice coefficient for active control servovalve = C^Wgy ~\f 2g c /p, 
0.00268 m 3 • sec' 1 • lb“^/^ (344.4 in 3 /sec / "^ lb / in®') 

f Coulomb friction between shock strut piston and cylinder, N (lb) 

F a vertical force exerted on shock strut by the runway surface, N (lb) 

Fj| impact phase limit force, N (lb) 

F lim limit force, N (lb) 

F g shock strut force, N (lb) 

F W g wing /gear interface force, N (lb) 
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SYMBOLS (Continued) 


g 


Sc 




*a 

Kf 

k FDGE 

Ksv 

Kt 

K x 

L 

M 


M L 

My 

PE t 


acceleration due to gravity, 9.80 m/sec 8 (386 • in/sec 8 ) 

gravitational acceleration constant, 1 kg • m • N' 1 • sec 8 
(12 slug • in • lb -1 • sec' 8 ) 

input signal to electronic compensation networks, A 

output signal from electronic compensation networks, or input signal to 
active control servovalve, A (+0.040 A maximum) 

amplifier gain in active control loop, 0. 000040 A/V 

position feedback gain in strut position control loop. 196. 9 V/m (5.0 V/in) 

fraction of total strut stroke assumed available when computing impact 
phase force from equation 37 

position gain of servovalve in active control loop, 0. 0635 m/A (2. 50 m/A) 
effective spring rate of tire, 245 000 • N/m (1400 lb/in) 
gain in strut position control loop, 1 . 0 m/m (1. 0 in/in) 
total life force, N (lb) 

mass of airplane per gear, 1456 kg (99.8 slugs) 

mass of upper portion of landing gear (cylinder plus orifice plate attachment, 
kg (slug) 

mass of lower portion of landing gear (piston plus tire), 32. 1 kg (2. 20 slugs) 
upper mass, 1456 kg (99.8 slugs) 

potential energy stored in tire due to compression, N-m (ft- lb) 
hydraulic supply pressure, 2.07 (10 7 ) N/m 8 (3000 lb/in 8 ) 
hydraulic return pressure, 0.0 N/m 8 (0.0 lb/in 3 ) 
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SYMBOLS (Continued) 


PROM programmable read only memory 

P x hydraulic pressure in shock strut piston, N/m 2 (lb/in 2 ) 

P s pneumatic pressure in shock strut cylinder, N/m 2 (lb/in 2 ) 

P 3 pressure in volume between walls of shock strut piston and cylinder, 

N/m 2 (lb/in 8 ) 

Q 0 flow rate through shock strut orifice from piston to cylinder, m 3 /sec 

(itf/sec) 

Q sv flow rate from active control servovalve to shock strut piston, linear 
mode, m 3 /sec (in 3 /sec) 

Q SVl flow rate through active control servovalve from supply pressure to the 
shock strut piston, m 3 /sec (in 3 /sec) 

Q SV2 flow rate through active control servovalve from shock strut piston to 
return pressure, m 3 /sec (in 3 /sec) 

RAM random access memory 

R s the slope of the limit force with respect to time during transition phase, 

137 900 N/sec (31 000 lb/sec) 

s LaPlace operator, sec' 1 

t time, sec 

V velocity 

V x hydraulic volume in shock strut piston and lines up to the active control 

servovalve, 0. 00426 m 3 (260. 0 in 3 ) for fully extended strut 

V s pneumatic volume, 0.000624 m 3 (38.1 in 3 ) for fully extended strut 

V a volume between shock strut piston and cylinder, 0.0 m 3 (0.0 in 3 ) for 

fully extended strut 


Wsv window width of orifices on third stage spool of active control servovalve, 
0. 0884 m (3.48 in) 
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SYMBOLS (Continued) 


displacement, m (in) 

commanded position of shock strut, 0. 1016 m (4.0 in) 
bulk modulus of hydraulic fluid, 6. 89 (10®) N/m a (1*10 lb/in ) 
ratio of specific healt of gas at constant pressure to that at constant 
volume, 1.1 

mass density of hydraulic fluid, 838 kg/m 3 (0.000941 slugs/in 3 ) 
time constant in strut position feedback loop, 0. 10 sec 
time constant in electronic compensation network, 0. 0281 sec 
time constant in electronic compensation network, 0.0141 sec 
time constant in electronic compensation network, 0.0010 sec 
time constant in electronic compensation network, 0.0001 sec 

corner frequency in active control servovalve transfer function, 

1263 sec -1 

natural frequency in active control servovalve transfer function, 

655. 5 sec -1 

natural frequency in electronic compensation network, 251.3 sec 1 
damping coefficient in active control servovalve transfer function, 0.436 
damping coefficient in electronic compensation network, 5. 10 
damping coefficient in electronic compensation network, 0, 10 


SYMBOLS (Continued) 


Subscripts: 

a lower mass of shock strut or axle 

g ground 

i initial conditions before impact 

im impact phase 

L lower mass 

max maximum value 

min minimum value 

r rollout phase 

s shock strut relative motion of the lower mass (piston) with respect to 

the upper mass (cylinder) 

sv servovalve 

tr transition phase 

U upper mass 

wg wing/gear interface 

Miscellaneous: 

d( ) indicates the differential of a variable 

) indicates difference or change in a variable 

• • • Ml 

( ),( ),( )dots indicate differentiation with respect to time 



ANALYSIS 


Nonlinear Math Model for Vertical Drops 

A simplified nonlinear dynamic model of an active control landing gear 
for vertical drops without aerodynamics and airplane dynamics was developed. 
The airplane mass is represented simply as a point mass located on top of 
the landing gear. This model proved to be of great value in this project for 
various reasons. First, it enhanced understanding of the basic equations of 
motion of the landing gear and the dynamic characteristics. This was very 
useful in the further development of an even more simplified linearized 
dynamic model, which was used extensively in the development of the elec- 
tronic compensation networks that are required for stability. The linear 
model and its application are discussed in later sections. Secondly, it often 
was desirable to study the dynamic operation of the system without the com- 
plicating influences of the airplane dynamics, and the simplified nonlinear 
model provided a good tool to do this. Also, by neglecting airplane dynamics, 
the size of the computer simulation is greatly reduced and runs can be made 
in less computational time. Finally, the simplified model was used as a 
check with the more general computer program supplied by NASA, which 
required many programming changes in its application to this study. 

Figure 1 shows a schematic of the simplified vertical drop case along 
with some of the variables used in modeling the system. The system basi- 
cally consists of (1) an upper mass M-^, which is the mass M (per gear) of 
the airplane plus the mass M^; of the upper portion of the landing gear (i. e. , 
the cylinder plus the orifice plate attachment), (2) a lower mass M L , which 
consists of the tire mass plus the mass of the lower portion of the landing 
gear (i. e. , the piston), (3) a three-way servovalve that ports fluid in the 
piston chamber to either supply or return pressure, and (4) an electronic 
controller. Note that the tire dynamics are simulated as a mass and linear 
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TOTAL 
UPPER MASS 


= M u - M + M c 



Figure 1. - Illustration of variables used in nonlinear 
simulation of simplified vertical drop case. 
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spring Kj. . The controller dynamics are presented later and will not be re- 
peated here. Suffice it to say that at each instant in time the controller 
processes the input signals to provide a current to the servovalve. 

Equations of Motion 

The equations of motion of the landing gear in the vertical direction will 
now be derived, with down being defined as the positive direction. The two 
independent masses and give rise to two degrees of freedom, X W g 
and X a (See Figure 1 ). The final results, however, will be expressed in 
terms of the displacement variables X w g and X g , where X g is the shock strut 
stroke, which ranges from zero for a fully extended strut to the maximum 
value for a fully contracted strut. The variable X s is equal to (X w g-X a j. 

Referring to Figure 2 , a force balance on the airplane mass gives 

MX wg « Mg - L-F S =-F wg ... (1) 

where F g is the "shock strut force" and F wg is the "wing gear interface 
force". 

Referring to Figure 3 , a force balance on the upper portion of the shock 
strut, which is assumed to be rigidly connected to the airplane mass, gives 

M c X w g = M c g + F s - P 2 A 2 + P2 (Ai-A 0 ) - Pi (Ai-A 0 ) + P3A3 + f 

or 

F s = M c X wg - M c g + (Pi-P2)(Ai-A 0 ) + P 2 A 2 - P3A3 + f ...( 2 ) 

Referring to Figure 4 , a force balance on the lower portion of the shock 
strut, which includes the wheel and tire mass, gives 
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L (LIFT) 


cn 
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Figure 3. - Free body diagram Figure 4. - Free body diagram 

of upper shock strut of lower shock strut 



• • • 


(3) 


ML,X a = M^g + P 2 (A2-A1 ) + Pi Ai - P3A3 - F a + f 

substituting (X W g-X g ) for X a 

M L (X wg - X S ) = M L g + ( Pj - P 2 ) Aj + P 2 A 2 - P3A3 -F a M 

where F a = ( (X wg - X s - X g ) for X wg - X s > X g 

I 0 for X wg - Xg ^ X g 

Note that the tire force F a is zero when the tire is not in contact with the 
runway surface. 

Equations were developed for: (1) the case where the piston is in motion 
relative to the cylinder, or V s 4 0; and (2) the case where the piston is not in 
motion relative to the cylinder due to Coulomb stiction between the two, 
or V = C. For the first case, the Coulomb friction f is simply the value of 
running friction that exists. But for the second case, the Coulomb friction is 
actually a stiction force that can vary over a r Singe of absolute values from 
zero up to the breakaway friction value. Transition from case 1 to case 2 
will occur whenever the relative velocity is zero or passes through zero 
throughout the digital integration process. Transition from case 2 to case 1 
will occur whenever the absolute value of the friction force f exceeds the 
known breakaway friction value. The following equations express these re- 
lationships. 

First, for the case when V a 4 0, equations 1 and 2 are combined to 
eliminate F s , 

M u X wg = M u g - L - (P 1 -P 2 )(A 1 -A 0 ) - P 2 A 2 + P3A3 + f ... (5) 
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Substituting equation 5 for X into equation 4 

wg 


M L 


M 


u 


[M u g - L - (P 1 -P 2 )(A 1 -A 0 ) - P 2 A 2 + P3A3 + f] - M L X S 


= M L g + (P 1 -P 2 )A 1 + P 2 A 2 - P3A3 - F a + f 


Solving for X g 

Mt Ml 

Ml Xs = [(Mu + M D g - L] + (1 + j^)[-M L g - (Pi-P 2 ) A l - P 2 A 2 + P3A3] 


m l 


Ml, 


+ M" (p l - P 2> A o + F a + (1 +^> f 


U 


M, 


u 


or 


M t X a 

L, S 


M L 

(1 )[FORCE~ f] 


•••( 6 ) 


Mt M y IVI t 

mT [ (Mu + Ml) g * L 1 + (1 + MT ) t‘ M L g ■ (P 1 " P 2> A 1 - p 2 a 2 + p 3 A 3) + TT“ ( p l - p 2> A o +F a 
FORCE = S — M u 


1 + 


M 

M„ 


L, 


Equations 5 and 6 constitute the two inertial equations of motion when V g f 0. 
Next, for the case when V g = 0, equations 4 and 5 are totaled to get 


(M u + M l ) X wg = (M u F M l ) g - L + (P t - P 2 ) A 0 - F a • • • (7) 


X s = 0 ...(8) 

Equations 7 and 8 are the two inertial equations of motion when V s = 0. 
However, an expression must be derived to tell when the two masses break 
away from one another; i. e. , when the break away stiction will be overcome 
causing relative motion between them. This required solving for the 
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variable friction force f. 


Substituting equation 5 for X 


wg 


into equation 4 


Mr 

-P [M u g - L - (Pi - P2KA1 - A 0 ) - P2A2 + P 3 a 3 + f l 
M u 

= M L g + (P x - P 2 ) A ! + P 2 A 2 - P 3 A 3 ‘ F a t f 


Solving for f, gives 


+ f = 


^ [(M u + M l ) g - L] + (1 + i^)| - M L g - (Pi - P 2 ) Ai - P 2 A 2 *• P3 a 3) 
— — — " 


m l 


If the value of f as calculated by equation 9 exceeds the breakaway stiction, 
then relative motion will result. 


Pressure, Flow, and Volume Relationship 

Referring to Figure 1, conservation of mass applied to the hydraulic 
fluid in volume Vj gives 


—£• Pi = A i X s - Q 0 + Qsvl - Qsv2 


...( 10 ) 
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where 


Qr 


Ao C d oV (Pi - P 2 ) for Pi > P 2 


'Aq *“do 


2 §c 


(P 2 - Pi) for P 2 > Pi 


V 1 = V^-AiXa 


The parameters Q gv ^ su^d Q are the servovalve flows from supply 
pressure to and from P^ to return pressure, respectively, and are a 
function of the spool displacement and valve geometry. These flows are cal- 
culated in a subroutine titled FLOZE2 in the nonlinear simulations. Note 
also that equation 10 considers the effect of fluid compressibility, which gives 
rise to the first derivative of pressure term. 

The pressure /volume relationship for the pneumatic volume V„ is 

a 


v 2[ ' 

p 2 = p 2 • < ) 

* V 2 


...dl) 


where 

v 2 = V 2i - (A 2 - Ai) X s - V cum + (V 3 - V 3 .) 

V CU m = ^ (Q q At) = the cumulative fluid flow through the 

orifice from Vi to V 2 

V 3 = V 3i + A 3 X s 
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The initial values of the pressure and volume variables referred to in the 
above equations are for the fully extended strut just prior to impact. 


Servovalve Dynamics 

The dynamic characteristics of the NASA-selected three-stage 
high response servovalve were supplied to HR by NASA. A mathematical 
representation of the servovalve dynamics was derived by the following pro- 
cedure. The small signal frequency response characteristics obtained from 
NASA are shown reproduced in Figure 5 for reference. Amplitude and phase 
angle data at several discrete frequencies over the range of interest were 
taken from these curves and input to an HR -developed computer program 
that solves for an approximate linear transfer function by the method of least 
squares. The resultant transfer function is 


X S v 

l 2 


K 


sv 


_f!_ + il§y s+ 1)( _1 + i) 

w8 sv W sv 


...( 12 ) 


where the values of the constants K gv> W gv * w c » and C gv are given in the 
SYMBOLS section and ig is input current and X gv is output displacement of 
the third stage spool. For the nonlinear simulation this transfer function is 
represented in differential equation form as follows: 
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i ••• , l 2 £sv 

oF - * x sv + *<£2 + w 

sv c sv 


0 ) 

sv C 


■) X sv + (tv- SV + ) X sv + X sv - K sv i2 

w sv c 


...(13) 


Saturations in the various stages are treated as follows. Third stage 
spool saturation obviously places a limit on X gv equal to ± X gy max. Second 
stage saturation is assumed to place a limit on d(X gv )/dt equal to ± 
(QSV2/ASV3), where QSV2 is the second stage saturation flow and ASV3 is 
the effective cross sectional area of the third stage spool connected to the 
lines from the second stage. First stage flapper saturation is assumed to place 
a limit on d 2 <X gv )/dt 2 equal to ± (QSV1* CQ 2 )/(ASV2-ASV3), where QSV1 is 
the first stage saturation flow, CQ 2 is the linearized flow gain of the second 
stage spool (flow rate per unit displacement of spool), ASV2 is the effective 
cross sectional area of the second stage spool connected to the flow lines 
from the first stage, and ASV3 is as already described. 

Linear Math Model for Vertical Drops 


This section describes the development of an approximate linear model of 
an active control landing gear for vertical drops, exclusive of aerodynamics 
and airplane dynamics. The linear model is a valuable tool since the effect 
of system modifications can be ascertained very rapidly. The study of system 
stability and the development of the electronic compensation networks, dis- 
cussed in the following section, relied very heavily on the linear model. The 
results were then supported by analysis using the complete nonlinear model. 

The equations for the actual landing gear and controller are inherently 
nonlinear, so some assumptions and simplifications must be made. The 
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r;« " — *“ — — - «- . 

«• «- - »• * u 

(2) The pneumatic pressure P 2 is assumed constant. 

(3) The pressure Pg is identical to the pressure P . 

(4) The lift and gravitational forces are constant. 2 

to thl 5 l The „° rOSS SeCtl0nal area ° f the PiSt0n 18 negli ® ible wi ‘ h aspect 
to the overall piston area. 

Piston 1 ^ ° f ‘ he ° rMCe “ " egligible wtth res ^ to the overall 
piston area. 

(7) The inertia of the upper part of the shock strut mass is negligible 
with respect to the airplane mass. 

(8) Servovalve nonlinearities are neglected (i & __ . , _ 

neglected (i.e., constant flow gain and 

pressure gam partial derivative values are assumed). 

The equations given in the previous section for the nonlinear math model 
can now be simplified as follows. Equations 5 and 6 become 


M n X 
U Wg 


M|jg - Li - PjAj 




“ L Xs “ M u l(M “ + M U g - L] + (1 + [- ML g - P lAll + K t (X wg - x s ) 


If we let the variables X 1 . 

wg* 

variables of X , X X 
wg* a* s’ 
variables about some mean 

written 


X a* X s* and p i represent the perturbation 
and P r (that is. they represent changes in those 
condition), then these two equations can be 
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- P'lAl 


M u X' wg 


M L X' S ■ < 1 + ^ )( - p, l A l )+ K t «' wg -X' s ) 


where the lift and gravitational terms have dropped out because they are 
assumed constant. Making use of the identity X Q = - X & , these two 

equations can be solved for X wg and X & . Dropping the primes for simplicity, 
the results in LaPLace transform notation are 


X wg 


AjPl 

M u s 2 


...(14) 


AjPl 

M L s 2 + K t • • *^ 15 ^ 


Also, note from equations 1 and 2 that 


F wg = F s « PiAi 


(16) 
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H: is instructive to note that equations 14 and 15 are the linearized equa- 
tions of motion that would result from a mass and piston driving into a sleeve 
connected to ground by a spring, with an orifice in the piston connecting the 
internal pressure Pj with a constant pressure P 2 on the other side. This is 
illustrated schematically in Figure 6, along with a servovalve for active 
control. A force balance on the masses M and M will result exactly in 
equations l4t and 15. 

A linearized expression for the pressure P^ will now be derived. 
Referring to Figure 6, conservation of mass applied to the hydraulic fluid in 
volume Vj gives 


dQsv “ d Qo -^i x c 


Yl p 
P 1 


...(17) 


The flow from the servovalve is 


*“SV x sv 


Q 


sv 


\fps - Pi 

0 


for X sv = 0 


^SV ^SV \ / P 1 - P R fo ^ x S v < 0 


Differentiating this, gives 


• sv 




g - Pj dX sv 


Csv Xsv 

2 VPs - Pi 


dPi 


for X sv ^ 0 


dQ 


sv 




C SV 1 / P 1 - P R dx sv + 


^SV x sv 


T^T 1 


dPi for X sv < 0 
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Letting the mean value of the pressure P equal (P + p )/2, and 

^ SR 

assuming small perturbations about this mean condition, the above relation- 
ships can be expressed as 


d Q S v “ CQ • dXgy - CP • dPj 


...(18) 


where 


CQ 


dQ sv I 

&X S y = C sv iA p S+Pr)/2 


CP 


SQsv 


- CQ/PG 


and where PG is the pressure gain of the spool in terms of pressure per 
spool displacement, at constant flow. In the linearized model the values of 
CQ and CP are assume constant. 

The flow through the fixed orifice is 



Differentiating this gives 


dQ 0 


K r 


2 V p l - p 2 


dP! a CP G dP x 


...(19) 
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However, note that at the mean condition Pj. equals P 2 , so the coefficient 
CP is undefined. In actuality, CP q is a function of excitation amplitude and 
frequency, ft is important to keep this fact in mind when applying the linear - 
ized model. The approach taken in this study was to vary CP q over a reason- 
able range to determine its influence on the dynamic response of the system, 
and to proceed accordingly . 

Substituting equations 18 and 19 into equation 17, solving for P^, and 
expressing it in LaPlace transform natation gives 


P 


1 


CQ • Xgy + Ajs Xg 

^ s + (CP 0 + CP) 
P 


( 20 ) 


The transfer function of the servovalve dynamics in terms of output dis - 
placement of the third stage spool for input current was derived in the pre- 
ceding section (equation 12) and is 


53V 

h 


K 


sv 


2 Csv , 
+ s + 1 




sv 


03 


sv 


s 


+ 1 


The transfer function for the electronic compensation network, which is 
derived in the next section, is 


y 2 

h 


/ S 2 + 2C2 k'is + u>iA I TjS + l \ 

\s 2 + 2d W 1S + “1 2 ) \ T 2 S + 1 ) 


\ r 4S+ 1 ) 


...( 21 ) 
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The relationship describing the force feedback and the position loop 
feedback is 


‘l - K a (Film - Fwg) - K a (^fj) (X c - K x X s > . ^ 

Equations 14, 15, 16, 20, 12, 21 and 22 make up the complete linarized 
model, and Figure 7 is a block diagram arrangement of these equations. The 
algebraic solution of transfer functions from this set of equations would be 
very tedious, so this practice was avoided. The following approach was 
taken. First, working from the block diagram, the system of equations is 
expressed in matrix form, where each element of the matrix is in general a 
polynomial in "s". A computer program is then used to solve for the co- 
efficients of the numerator and denominator of the desired transfer function 
using Cramer's Rule. Frequency response characteristics are then calculated 
from the transfer function. Also, the roots of the numerator and denominator 
polynominals are easily determined using standard computer techniques. The 
linear model gives very rapid results, and thus is ideal for evaluating the 
effect of varying different parameters and for developing electronic compen- 
sation networks for stability. 

The validity of the linear model was evaluated by comparing the open 
loop frequency response with results obtained from the nonlinear vertical drop 
model described in the preceding section. The loop was opened at the point of 
wing gear force (F ) feedback, and position feedback was not included. 

The conditions were for zero lift with the shock strut at a mean stroke equal to 
that required to balance the weight of the airplane, which was 67 percent of 
maximum stroke for this case. Also, no electronic compensation was used. 
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Figure 7. - Block diagram of linear math model 










The comparison of results is shown in Figure 8. Input is commanded limit 
force and output is wing gear force. Linear model results are shown for both 
a low and a high value of the fixed orifice parameter CP q . The results for the 
low value of CP Q match the nonlinear results closest at low frequencies while 
the results for the high value of CP q match closest at high frequencies. The 
basic trends are predicted reasonably well and it is concluded that the linear 
model is a valid tool for performing initial design and stability studies.. 


Energy Considerations 


The primary function of the landing gear control system is to control the 
wing gear force at some commanded level. A secondary but necessary 
function of the controller is to set this commanded level such that the strut 
will not fully collapse, precluding structural damage. In order to do this, 
the controller must be able to monitor the energies present in the system during 
impact, and then set the limit force accordingly so that the work potential of 
the shock strut will be sufficient to absorb these energies. In order to study 
the feasibility of the controller performing this function, it is worthwhile to 
undergo an analytical development of the various energies present in the 
system and relate them to the work performed by the shock strut — the 
subject of this section. The equations will be developed from the equations of 
motion of the simplified vertical drop case, presented earlier; however, the 
results are almost directly applicable to the more general cases. 

Equations 3 and 5 are the inertial equations of motion corresponding to 
the two independent masses M_ and M , These equations are rewritten in a 

U 1-4 

form where each term represents the differential of an energy or work. This 

is done by multiplying equation 5 by dX and equation 3 by dX , resulting in 

wg a, 

M u V wg dV wg = [M u g - L - (p l - P 2>( A 1 - A o> - P 2 A 2 + P 3 A 3 Tf] dX wg 

... (23) 
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re 8. - Comparison of linear model and nonlinear model results 
open loop, no compensation, no position feedback. 










M L V a dV a = [M L g + P 2 ( A 2 ~ A l> + p lAi - P3A3 - F a + f] dX a ' 

...( 24 ) 

Note that if these equations are integrated from some initial condition to a 
final condition they express changes in energy that occur over that time 
interval. For example, the integral of the terms on the left side of the equa- 
tions express the total change in kinetic energy of each mass. It is desired 
now to write an equation expressing the total change in kinetic energy of the 
system, which is assumed to consist entirely of the two masses M u and M 
as a function of the various other energy and work terms. This is done 
by adding equations 23 and 24, resulting in 


1 

✓Si 


2 


M u v wg dV wg + M L v a dV a + [~M u g + L -(Pi-P 2>A 0 ] dX. 


wg 


9 


10 


+ f _M LS + F al dx a + K p i ~ p 2^ A 1 + P 2 A 2“ P 3 A 3 ± f) dX 


= 0 


. . • ( 25 ) 

Equation 25 can be thought of as a statement of conservation of energy, 
since it essentially says that the sums of all changes in energy and work done 
in the system equals zero. Some of the terms in equation 25 are easily 
recognizable as to the type of energy they represent. For example, terms 
1 and 2 represent changes in kinetic energy. Terms 3 and 6 represent 
changes in potential energy of the upper and lower masses due to height 
changes. Term 7 represents changes in potential energy stored in the tire 
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spring. Term 10 represents energy dissipation due to Coulomb friction. 
Energy dissipation due to flow across the shock strut orifice also occurs, 
although this is not directly identifiable in equation 25 because of the com- 
plicating effects of fluid compressibility and flow to and from the active - 
control servovalve. 

Equation 25 must be simplified to the point where it becomes a useful 
tool to the controller for monitoring the energies present in the system in 
terms of the available input signals; e. g. , accelerometers, transducers, etc 
First, term 5 was found to be negligible and can be immediately eliminated. 
Next, the combination of terms in the third bracket (terms 8, 9 and 10) are 
approximately equivalent to the shock strut force F g , This can be seen by 
referring back to equation 2, neglecting the inertia and gravity forces on the 
shock strut portion of the upper mass, and neglecting the orifice area with 
respect to the piston area. With these simplifications, equation 25 is re- 
written as 


MuV w gdV we + M L V a dV a + (-M ug + L) dX wg + (-M L g+ F a ) dX a + F g dX s = 0 

...(26) 

Note that the quantity F g dX g represents the differential of the work performed 
by the shock strut. Integration of this equation over a given interval will thus 
yield an expression showing the net work performed by the shock strut as a 
function of the various other energy changes. The task of the controller is to 
monitor these energies and make sure that the remaining work potential of 

the shock strut is sufficient to dissipate them. 

At this point, it is worthwhile summarizing the instrumentation that is 

available to the controller as originally proposed. 

• • 

Cl) Wing/gear acceleromters - measures X 

•9 & 

(2) Hub accelerometer - measures X a 

(3) LVDT on shock struct - measures X g 

(4) Upper cylinder pressure transducer - measures P 2 (pneumatic) 

(5) Lower piston pressure transducer - measures P 1 (hydraulic). 
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With this instrumentation, it is not possible to apply equation 26 to the task 
at hand, because there is no way to monitor the lift L or the shock strut 
force F g (note that by knowing one of these, the other may be obtained 
equation 1), It is therefore necessary to make further simplifications. 
Assuming that the total airplane weight equals the lift throughout impact, and 
neglecting the M^g term, then equation 26 becomes 

MuV wg dV wg + M L V a dV a + F a dX a + F g dXg = 0 

...(27) 

With the weight -equals -lift assumption, all of the parameters in equation 27 
may be obtained. From equation 1, 

F s =_M uX wg ...(28) 

It is possible to solve for F g by first noting that equation 3 may be expressed 
as 

Mj_X a *=* -F a + F s 
Substituting equation 28 for F g into equation 29, 

F a ~ - Mu X - M l X a 

...(30) 

The individual terms in equation 27 were integrated with respect to time 
throughout a typical vertical drop simulation using the simplified nonlinear 
vertical drop model described in a previous section to determine if the 
equation is indeed a valid representation of all the major energy changes that 
take place. The lift was set equal to the weight throughout the simulation, so 


...(29) 

f 
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this test case will say nothing regarding the validity of the "weight -equals - 
lift" simplification. The validity of that will be demonstrated later using 
actual land and roll simulations. The integrated variables are shown plotted 
in Figure 9. The sum of all the integrated variables is also shown, repre- 
senting the total energy plus work done by the system. This should remain 
constant from the conservation of energy principle, ft is indeed relatively 
constant, thus supporting the validity of equation 27, at least for vertical drop 
cases where weight equals lift. Note that the kinetic energy of the lower mass 
is negligible. Also, the potential energy stored in the tire, although not 
completely negligible, is a minor part of the total energy of the system. 

This raises the question whether it is worth the added complexity of including 
this term in the energy computations of the controller. The approach taken 
here is to neglect it and incorporate any errors into an empirical correction 
factor, if necessary. As a result, the only terms considered in the energy 
computations of the controller are the kinetic energy of the upper mass and 
the work of the’ shock strut. 

The next step is to derive relationships that the controller can apply in 
setting the commanded limit force to a level such that bottoming of the strut 
will be prevented. There are various approximations or simplifications that 
can be made, some of which have already been mentioned in the preceding 
paragraph, which will introduce errors in the final result. In the interest of 
completeness and for reference, the more general derivation will be present- 
ed first, followed by the simplifications. The procedure is to integrate 
equation 27 from a time (controller initiation) to the end of transition; 
i. e. , the start of rollout. The following assumptions are made. 

(1) The actual wing/gear force is identical to the commanded limit force 

throughout the impact and transition phases. 

(2) The velocity of the masses is zero at the end of transition (this 
also is the criterion used in calculating the wing/gear velocity at which 
transition begins). 
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Figure 9. - Various energy terms versus time, from the 
simplified nonlinear vertical drop simulation. 














(3) The commanded limit force is constant throughout the impact phase 
and linearly decreases during transition to zero limit force at the beginning 
of rollout. The slope in the transition phase is a known quantity. 

A sketch of the variation of wing/gear force versus time for the above 
conditions is shown in Figure 10. Integration of equation 27 from t^ to the 
start of transition t^ gives 


1 /2 M u (V 2 - V 2 ) f 1 12 M l (V? - V 2 l ) 


tr 


tr 


L tr 

7 


F a dX a + F U AX S . = 0 


’im 


C 1 


Integration from the start of transition t tr to the end of transition t gives 

t. 


1/2 M u (y£ r - V^ r ) + 1/2 m l (v£ r - v[ tr ) + J 


F a dX a + ^- AX Str = 0 


tr 


The quantities AX s * m and represent the shock strut stroke used in the 

impact and transition phases, respectively. The sum of these should be 
equal to or less than the total remaining available stroke, to prevent bottom- 
ing. Solving for AX S . and AX S+ in these equations and adding them together 

l °xr 

gives t tr 

1/2 M U (V 2 -V 2 )+1/2M l V 2 - / F a dX a 

u u tr L tr J 

h 


AX 


s 


F 


li 


1/2 M U V 2 Utr + 1/2 M l V 2 l 


tr 


/; 

Hr 


F a dX a ...(31) 


F u /2 


Referring back to Figure 9, it is shown that any potential energy changes 
in the tire during the impact phase may be neglected. This is true because F a 
and X a remain relatively constant throughout this phase. Also, since the final 
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ACTIVE ACTIVE / TRANSITION 


UJ!l =i 'aoyod ±iiAin 
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Figure 10. - Illustration of the variation of limit force versus time 
and the various phases of operation following impact. 


rollout limit force is zero, F at that point will be approximately zero; and 

St 

hence, roughly stated, that during transition the potential energy of the tire 

wfii go from its initial impact phase state to a zero energy state. This fact 

can be seen in Figure 9. Thus, in equation 31, the integral of F & dX a over the 

interval t, to t. is approximately zero and the integral of F dX over the 
1 tr a a 

interval t. to t is approximately equal to the negative of the potential energy 
tr r 

stored in the tire at time tj, the start of the integration interval. Signifying 
this energy level by the symbol PE^., and neglecting the kinetic energy of the 
lower mass, equation 31 reduces to 


F u £X S = 1/2 M u V 2 u + 1/2 M u (V Utr ) 2 + 2 • PE t 

...(32) 


The upper mass velocity at the start of transition can be determined by equat - 
ing the impulse to the change in momentum for the transition phase, resulting 

in 

F 2 

V = — 

u tr 2 • M u • R s 

u s ...(33) 

where R is the slope of the limit force with respect to time during transition, 
s 

The application of equation 32 to the controller could be as follows: 
when the gear first impacts the runway surface the controller is in a passive 
mode; that is, no command is being initiated to the servovalve. Immediately 
after impact, at discrete intervals, the controller will compute the strut work 
(F^AXg) required to absorb the energy using equations 32 and 33. In equa- 
tion 33, the transition velocity is calculated using the current value of wing/ 
gear interface force as the impact limit force (i.e., F ^ is set equal to F ). 
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This calculated value of required strut work is continually compared to a 
signal representing the actual strut work currently available if the controller 
were to go into active control at that instant. This actual strut work is simply 
the current value of the wing gear force times the remaining shock strut 
stroke. When comparison of these signals indicates that the available shock 
strut work exceeds the value calculated by equation 32, then active control to 
the servovalve is initiated and the commanded impact phase limit force is set 
equal to that current value of wing gear force. The velocity of transition is 
also then set to the value computed from equation 33. 

It may be desirable to further monitor the energy situation through the 
remainder of the impact phase to ensure that the strut will not bottom due to 
additional energy inputs. If so, equations 32 and 33 could be used to compute 

from the known remaining strut stroke (AX g ), and if the computed F^ 
ever exceeded the original F^ value, the controller would update F^. to the 
newly -computed value; however, this would require the solution of a quartic 
equation. 

A simplification in equation 32 can be made by neglecting the transition 
region. In other words, assume that in Figure 10 the constant impact limit 
force is continued throughout the transition region, and at the end of transi- 
tion it drops abruptly to the rollout limit force of zero. It should also be 
stressed that it is not being proposed that the actual limit force commanded 
by the controller follow this pattern; this is only for the purpose of simplify- 
ing the energy relationships. Integration of equation 27 will then result in 


1/2 M u (V 2 Ur - V 2 u ) + 1 /2 Ml (V 2 - V 2 ) + / F a dX a + F u AX S = 0 


r L 


/ 


...(34) 


37 



Again, the integral of F & dX a from * j_ to t r is approximately equal to the 
negative of the potential energy stored in the tire at time tj. Neglecting the 
kinetic energy of the lower mass, equation 34 becomes 


F U AX S = 1 12 M u V 2 u + PE t 

...(35) 

with equation 33 still applying for the velocity of transition. The use of 
equations 35 and 33 in the controller would be similar to that already described 
for equations 32 and 33. 

A still further simplification, in addition to those already made, is to 
neglect the potential energy in the tire. Equation 35 then reduces to 


F u AX S = 1 12 M u V 2 


...(36) 


with equation 33 still applying for the velocity of transition. This last 
approach is the one taken in the current design, and is the one programmed 
into the computer simulations. 

To extend the use of equation 36 to actual land and roll cases, AX cos $ 

s 

is substituted in place of AX , where 6 is the angle the shock strut axis 
makes with the vertical (gravitational) direction. The angle $ is equal to the 
pitch angle of the airplane fuselage plus the angle of the wing with respect to 
the airplane fuselage. Also, the velocity and force in equation 36 are now the 
components in the vertical, or gravitational direction. The energy relation- 
ship used for land and roll cases is then 

F u Ax s cos0 = 1/2 M u V 2 u ...(37) 
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For land and roll cases, the question of the validity of the "weight-equals - 
lift" simplication is raised. Whether or not this is valid is demonstrated 
by computer simulation of the actual landing roll cases provided by NASA. 

The simulation results are presented in the ANALYTICAL RESULTS section. 

Of interest to the present discussion, however, is the maximum resultant 
stroke that occurs in each case when the limit force is set according to 
equation 37 and the aforementioned procedure. Table I, presented in a later 
section, summarizes this information. As can be seen, the method is reason- 
ably successful in setting the limit force so that a maximum amount of strut 
stroke is used in absorbing the impact energy. 

Stability Analysis 

The linear model discussed earlier was used to evaluate the frequency re- 
sponse and stability of the system. System parameters used were basically those 
given under SYMBOLS. The linearized orifice coefficient CP Q was given a value 
of 7. 13 (10 -7 ) m^/sec/kPa (0. 3 in^/sec/lb/in^) for frequencies less than 14 Hz 
and a value of 4.75 (10”®) m^/sec/kPa (2.0 in^/ sec/lb/ in^) for frequencies 
greater than 14 Hz, providing a reasonable match with frequency responses 
generated with the nonlinear model. 

Based on results from the nonlinear vertical drop model, it was determined 
that a force loop bandwidth of 10 Hz would be adequate to achieve desired per- 
formance. Using the uncompensated linear model shown in Figure 11, the 
Nyquist plot of Figure 12 was generated. The uncompensated system was 
unstable at a frequency (140 Hz) related to the servovalve dynamics and could 
not be stabilized by simply lowering loop gain without violating the 10-Hz 
bandwidth goal. It was evident that signal shaping was required. 

Analysis of various forms of forward path electronic compensation 
ultimately yielded the filter configuration shown in Figure 13. It consists of 
a second-order notch and two lead lags. The resulting Nyquist plot is shown 
in Figure 14; and, it can be seen that stability margins are adequate. 
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ire 12. - Nyquist plot of R^g/F^^ w ith force loop opened, 
,no position feedback, no compensation. 
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Figure 13. - Electronic compensation network and its characteristics 









































The strut position feedback was next incorporated into the models. Its 
purpose is to ensure the return of the strut to its neutral position after land- 
ing impact. With no specific requirements on the time to return to neutral, 
it was assumed that a relatively gradual (approximately 10-sec) return would 
be acceptable. This assumption allowed a low bandwidth position loop which 
is essentially decoupled from the force loop. The resulting benefit is that 
the position loop can be continually active at all times during impact, transi- 
tion, and rollout without degrading the performance of the force loop. A 
first-order lag in the position loop forward path yielded the desired results. 

The block diagram for*the complete linearized system ‘with position and force 
loops and compensation is shown in Figure 7. 

Figure 15 shows the Nyquist plot of the F w g/F^j^ m response for the 
complete system with the force loop opened and the position loop closed. 

Figure 16 shows the Nyquist plot of the X s /X c response with the position loop 
open and the force loop closed. Adequate stability margins are exhibited in 
these plots. 

Figures 17 and 18 show the frequency response characteristics of the 
complete linearized system with the force and position loops closed, for 
Fwg/^lim anc ^ X s /X c responses, respectively. Note that the bandwidth of 
the Fwg/Fiim response is about one hundred times greater than that for the 
X s /X c response. 

The sensitivity of the system to hydraulic fluid compressibility was 
investigated because the effective bulk modulus of the air/oil mixture within 
the strut is relatively unknown. Throughout this study, a nominal bulk 
modulus of 6. 89 (10 5 ) kPa (100 000 lb/ in 2 ) was used based on past experience 
with hydraulic systems. Nyquist plots indicating system stability with bulk 
modulus values of 3. 45 (10 5 ) kPa (50 000 lb/in 2 ) and 1. 38 (10 6 ) kPa (200 000 lb/ in 2 ) 
are shown in Figures 19 and 20. A linearized orifice coefficient of 
7. 13 ( 10“ 7 ) m 3 /sec/kPa (0. 3 in 3 / sec/ lb/ in 2 ) was used at all frequencies for 
these cases. A high bulk modulus is destabilizing in the 200-Hz range and a 


44 



























Figure 18. - Frequency response of X s /X c with force and position 
loops closed, with compensation. 












Figure 19. Nyquist plot of F w g/F^j m with force loop opened 
position loop closed, with compensation. 











low bulk modulus is destabilizing in the 100 -Hz range. While these results 
show some sensitivity to bulk modulus, they are thought to be conservative 
in light of the model being used. It is expected that the bulk modulus in the 
actual strut will lie somewhere between 3. 45 (10^) and 1* 38 (10 6 ) m^/sec/kPa 

(50 000 and 200 000 lb/in 2 ), and that the system, as configured, will exhibit 
the desired performance. 

ANALYTICAL RESULTS 

In this section, transient response results obtained from digital computer 
simulations are presented for: (1) vertical drop cases using the simplified 
nonlinear vertical drop model presented earlier, and (2) actual land and roll 
cases provided by NASA using a NASA -supplied computer program which was 
modified by HR for incorporation into this study. The purpose of these simu- 
lations is to demonstrate that the landing gear controller is developed and 
described in the previous sections performs the desired functions satis- 
factorily. 

The NASA computer program includes detailed simulation of the aero- 
dynamic and aircraft dynamics in two dimensions, as opposed to the one 
dimension for the simplified vertical drop model. However, the NASA pro- 
gram as supplied to HR. did not contain control system and servovalve dy - 
namics and did not include the effect of fluid compressibility. Therefore, it 
was modified to include all these items. 

Table I lists aH the cases which were simulated and reported in this 
section, along with a brief description of each. Cases will be referred to 
by the case number shown in the left column. For each case, two simulations 
were run one using an active control landing gear and one using a passive 
gear — so that the performance of the active gear in reducing wing forces 
could be evaluated. 
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SUMMARY OF ACTIVE CONTROL LANDING CASES 
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* K =0. 9 All other Cases K „„ 

FDGE FDGE 



Vertical Drop Cases 


Figure 21 shows the commanded and actual wing gear forces occurring 
during the initial impact and transition phases for Case A. The total lift 
force is set equal to the total weight throughout this transient; and, as a con- 
sequence, the landing gear eventually lifts off the ground (or rebounds) after 
0. 754 seconds. At that time, the gear is fully extended and the system has 
some finite upward velocity which it continues to maintain as long as the lift 
stays equal to the weight. Figure 21 also shows the forces for the passive 
gear case. A 28% reduction in peak wing/gear force is obtained when using 
active control. 

Figure 22 shows the forces for Case B. This is the same as Case A, 
except that starting from the instant of impact the lift force is linearly reduced 
to one half of its initial value. This reduction occurs over a time interval of 
0.4 seconds. The gear does not leave the ground for this case. A 23% re- 
duction in peak wing gear force from the passive gear case is obtained. The 
response of the strut stroke is shown in Figure 23. As the impact transients 
subside, the strut approaches the value commanded by the strut positioning 
loop. 


Land and Roll Cases 

Figure 24 shows the commanded and actual wing/gear forces occurring 
during the initial impact and transition phases for Case 1, which uses a 3 -Hz, 
0.0254 m (1-inch) amplitude sinusoidal runway for the first 3.083 seconds, then 
a flat runway. For this case, an 11% reduction in peak wing/ gear force is ob- 
tained when using active control. Figure 25 shows the pressures in the strut 
as a function of time. Figure 26 shows the strut stroke transient. At a time 
of 1.86 seconds, the strut hits the stops in the fully extended position and 
remains there until 2.01 seconds. Figures 27 and 28 show the wing gear force 
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Figure 21. - Wing/ gear force transient, case A, 













Figure 22. - Wing/ gear force transient, case B. 
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Figure 23. - Strut position transient, case 




Figure 24. - Wing/gear force transient. 
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Figure 25. - Strut pressures, case 





Figure 26. - Strut position transient 






















and strut stroke transients for Case 2, which is the same as Case 1 except it 
uses a flat runway and has no wheel braking, while Figures 29 and 30 show 
the wing/gear force transients for Cases 3 and 4. 

The flow rate through the servovalve of the active control landing gear 
during landing impacts is also of interest. Figures 31 through 36 show the 
time histories of the servovalve flow rates for each of the landings simulated 
herein. Positive flow is flow from supply to the piston chamber and negative 
flow is flow from the piston chamber to return. 

Table I summarizes the amount of wing gear force reduction that was 
obtained from using active control for each of the cases reported in this 
section. From inspection of the force transients it appears that the reduc- 
tion for some cases could be improved by increasing the bandwidth of the 
controller, which is presently about 10 Hz. This would reduce the magnitude 
of the initial overshoot. However, increasing the bandwidth would likely 
require more complex electronic compensation than that presently being 
used. In the absence of specific requirements, the results presented herein 
are judged satisfactory. The analyses indicate that the basic concept is 
workable. 
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Figure 29. - Wing/ gear force transient, case 3 
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Figure 32. - Servovalve flow rate, case B 






Figure 33. - Servovalve flow rate, case 










Figure 34. - Servovalve flow rate, case 2 

























Figure 36. - Servovalve flow rate, case 4 











CONTROLLER FUNCTIONAL DESCRIPTION 


Associated Equipment 

The following NASA- supplied equipment is used in conjunction with the 
electronic controller. 

(1) Accelerometer - range ± 3 g's, measures wing/ gear interface 
force. 

(2) Accelerometer - range ± 25 g's, measures hub acceleration. 

(3) Potentiometer - range 0-0. 305 m (0-12 in.), measures strut 
deflection. 

(4) Pressure transducer hydraulic - range 0-3.45 (10^) kPa 
( 0-5000 lb/ in ), measures the gear hydraulic pressure. 

( 5 ) Pressure transducer, pneumatic - range 0-1.38 (10^) kPa 
(0-2000 lb/ in ), measures the gear pneumatic pressure. 

(6) Servovalve - max. flow 0. 757 m^/min @ 6895 kPa 
(200 gal/min. @ 1000 lb/ in 2 ). 

(7) Servovalve Electronic controller - provides the signal to drive 
the servovalve. 

(8) Hydraulic power unit. 

(9) Modified landing gear. 

Interface 


The interface of input/ output signals is shown schematically in Figure , 
37. Primary inputs and outputs are those required by the controller to 
perform its function. Secondary inputs and outputs are provided for testing, 
status indication and parameter monitoring, but play no part in the basic 
function of the controller. The primary inputs are shown in Table II, 
secondary inputs in Table III and the outputs are in Table TV. 
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Figure 37. - Input/Output Signal Schematic. 




TABLE IL - PRIMARY INPUT SIGNAL SPECIFICATIONS 
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TABLE IH. - SECONDARY INPUT SPECIFICATIONS 
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TABLE IV. - OUTPUT SPECIFICATIONS 




Performance Characteristics 


The controller monitors the sensor data; and then, by performing 
computations of energy relationships and signal conditioning, provides an 
output signal which is the input command to the servovalve. 

The controller has three basic functions which are: (1) operating-mode 
determination, (2) limit-force command determination, and (3) control- 
law implementation. 


Ope ratin g- Mode Determination 

When the controller is enabled, it automatically determines the 
operating mode — landing or takeoff. 

Landing Mode 

The controller selects the landing mode if the controller -enable signal 
has been received and the strut position signal indicates full extension. The 
landing mode is divided into several phases, each imposing a different 
functional command on the controller. These are: (1) pre -touchdown, 

(2) active control initiation, (3) transition, and (4) rollout. 

Pre -touchdown. - During the pre-touchdown phase, the controller pro- 
vides a bias signal to the servovalve to maintain the strut hydraulic pressure 
equal to the design charging pressure. This is accomplished by a pressure 
control loop in which the hydraulic pressure is the feedback signal. It also 
receives a signal from an external source which is representative of the 
aircraft sink rate prior to touchdown. In addition, it monitors the strut 
deflection. The servoloop is not enabled and the gear remains in a passive 
state. 
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Active Control Initiat i on. - Active control is initiated when the energy- 
relationships indicate that the work potential of the strut exceeds the kinetic 
energy of the aircraft. Upon such occurrence, the controller causes the 
following to occur. 

(1) The servoloop is enabled. 

(2) An output is generated to the servovalve controller proportional 
to the force error as modified by the control laws. 

(3) Energy computations are discontinued. 

(4) A constant limit force is maintained. 

(5) The strut pressure loop is opened. 

(6) The servovalve bias is removed. 

(7) The transition velocity is computed and continuously compared 
to the actual wing gear velocity to determine the start of transition. 

Transition. - Transition to the rollout phase commences when the 
wing gear velocity equals the transition velocity. During transition, the 
controller linearly decreases the limit force command to zero and maintains 
active control. 

Rollout. - The rollout phase commences when the limit force command 
during transition reaches zero. Active control is maintained and the 
controller remains in this mode until receipt of a reset signal (push button 
on the front panel). 


Takeoff Mode 

The controller selects the takeoff mode when; (1) a controller-enable 
signal has been received, and (2) the strut potentiometer signal indicates 
that the strut is atO. 005 m (0. 2 in) less-than-full extension. At takeoff, the 
controller status is the same as in the rollout phase of the landing mode. 
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Limit -Force Command Determination 


During the period prior to control initiation, the limit force is zero. 
(Since the servoloop is disabled during this time, the limit-force command 
has no effect). During the active control phase, the limit-force command 
is equal to the wing gear force which was present at the instant when active 
control was initiated. As indicated above, the limit-force command is 
linearly decreased to zero during the transition phase and maintained at 
zero throughout the rollout phase. 

Control Law Implementation 

The controller implements the control laws as shown in Figure 38. The 
transfer functions are listed in Table V. 

Controller Design 
Functional Sections 

The controller consists of a digital section, an analog section and a 
power supply section as shown in Figure 39. 

Digital Section . - The digital section includes a general-purpose 
stored-program digital computer (microprocessor) and its associated input/ 
output devices as shown in Figure 40. The digital section accepts three 
analog and three discrete inputs and provides one analog and two discrete 
outputs. The analog inputs are applied through an analog multiplexer and 
A/D converter. The multiplexer is controlled by the computer software. 
The operating range of all inputs is 0 to +10 V. 

The controller reset input is generated by a push button switch on the 
front panel. This input is normally in a logic "0 " state. Depression of the 
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WING/GEAR ACCELEROMETER 
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Figure 38. - Control Law Functional Schematic 




TABLE V. - CONTROL LAW TRANSFER FUNCTIONS 
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Figure 39. - Controller Functional Schematic 











button changes the state to ”1" which causes the controller to reset the 
computer. (The computer goes to step 0000 of memory). 

The controller is enabled by a logic level input. Provision is made so 
that this input can be generated by a panel-mounted push button switch or a 
microswitch mounted on the gear. Both switches are connected in series. 
Therefore, if the controller is to be enabled by the panel switch only, then 
the microswitch input should be short-circuited. In any event, this signal 
is normally in a "0" state, and upon actuation the state becomes "1" which 
causes the computer to enable the controller. 

Initial sink rate is an additional analog input to the multiplexer. This 
input is generated by an external source which is a velocity sensor or a 
d c signal which simulates the sink rate. 

The outputs from the controller are the limit-force command (analog) 

signal, servoloop-enable (discrete) signal, and integrator -enable (discrete) 
signal. 

The limit-force command is an output from the digital section by means 
of a digital- to - analog converter. The servoloop-enable signal is a logic 
level signal. A logic level "1” causes a switch to close and completes the 
path between the limit-force command and the output to the servovalve con- 
troller. The integrator-enable signal is also a logic-level signal. A level 
"0" causes a switch which is connected across the integrator, to open, 
thus allowing integration to occur. The integrator is used to generate wing/ 
gear velocity from wing/gear acceleration and the input sink rate 

The computer software is illustrated in flow diagram form in Figure 41 
and a detailed listing is provided as Appendix a. 

A nalog Section . - The function of the analog section is shown schemati- 
cally in Figure 42. A detailed schematic drawing, wiring diagrams and 
assembly drawings, as well as a complete functional description can be 
found in Appendix B. The analog section accepts the primary and secondary 
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Figure 41. - Software Flow Chart 
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Figure 42. - Analog Section Schematic 










inputs, implements the required control laws and provides the output signal 
to the servovalve controller. It accepts signals from, and supplies signals 
to, the digital section. 

Provisions are included to apply test signals to the controller by means 
of jacks on the front panel. These are shown in Figure 43. Buffered outputs 
are also provided on the front panel so that signals of interest can be moni- 
tQred or recorded. These are shown in Figure 44. 

The inputs to the analog section from the digital section are the limit- , 
force command, servoloop- enable and integrator— enable signals, while the 
outputs from the analog section to the digital section are the wing/ gear 
velocity, wing/ gear acceleration and strut position signals. As previously 
stated, the wing/ gear velocity is derived by integrating the wing/ gear 
acceleration. 

In addition, the analog section accepts inputs from the following 
sensors. 

(1) Wing/ gear accelerometer. 

(2) Hub accelerometer. 

(3) Strut position potentiometer. 

(4) Strut hydraulic pressure transducer. 

(5) Strut pneumatic pressure transducer. 

These sig nal s are applied by means of connectors on the rear panel. 

Several other signals are applied to the analog section by means of 
controls on the front panel. These are: 

(1) Servovalve bias command - to determine the static 
operating level of the servovalve. 

(2) Strut position command - to determine static strut extension. 
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TEST 

Test Fixture 


The test fixture used for the ACLG is a pivoted, variabie -inertia beam 
originally designed to test the S-1C main engine vectoring servoactuators. 

The fixture is shown in Figures 45 and 46. Configured for the ACLG test 
program, the beam, added weights, landing gear, associated mounting com- 
ponents and miscellaneous hardware, represents a total calculated rotational 
inertia of 1426 kg - m2 (1050 slug-f t 2). The landing gear , g rigidly attached at 

a distance of 1 . 02 m (3. 33 ft) from the beam pivot which produces a vertical 
translational inertia equivalent to a mass of 1388 kg (3060 lb>. This value is 
about 5% below the targeted simulation value of 1458 kg (3215 lb) in order to 

allow for the weight of components which were added (hoses, fittings, instru- 
mentation, etc.). 

Motion is imparted to the beam by means of a hydraulic cylinder attached 
(through a load cell) to the beam at a distance of 0. 305 m (1 ft) from the beam 
pivot and on the opposite side of the pivot from the landing gear. The cylinder 
is controlled by a 0. 00189-m3 /s6c (30-gal/min) electrohydraulic servovalve 

and its associated electronics. An LVDT is used for position indication and 
position loop closure. 

Since the gear is rigidly mounted to the beam, it traverses the same angle 
as the beam. The angle of the gear at touchdown is about 4° forward of vertical. 
After touchdown, it passes through vertical and reaches an extreme of approx- 
imately 11° aft of vertical, depending on the test condition. These angles were 
chosen on the basis of gear geometry and the gear forces and deflections pre- 
dicted by the analytical simulation. They represent an optimum mounting con- 
figuration, the basis of which is to make the strut frictional forces the same 

percentage of the total axial force on the gear at both touchdown and maximum 
force condition. 
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Figure 46. - S-1C Test Fixture, Landing 
Gear Installed (View 2) 8 



The loading system operates in three modes as follows: 

(1) Pre-test gear positioning: The cylinder is in a closed position loop, 
using the LVDT for feedback to ’’cock" the beam and gear to a pre-test position 

as determined by a command potentiometer. 

(2) Sink speed generation: The LVDT is disconnected so that the cylinder 
is an open- loop configuration. A preset velocity potentiometer then provides 
a signal to the servovalve, thereby commanding a flow rate from the cylinder 

corresponding to the required sink speed. 

(3) Programmed force: At touchdown, a microswitch is activated which 

connects the load cell and places the cylinder in a closed force loop to control 
the force during the interval from touchdown to the end of the test. The force 
program is provided by a ramp generator which begins with a 4315 N (970 lb) 
upward force at the gear [14 370 N (3230 lb) of tension at the load cylinder] which 
cancels the beam unbalance of 4315 N (970 lb) and simulates the condition of lift = 
weight. After 0.4 sec, the force ramps linearly to provide a downward force 
on the gear of 2180 N (490 lb) [7250 N (1630 lb) of compression at the load cylinder 
which adds to the beam unbalance at this point of 4982 N (1120 lb) (due to the fact 
that the beam c.g. is above its pivot point) to produce a total downward force of 
7181 N (1610 lb), which is equivalent to half of the aircraft weight per gear. The 
force is then held constant at this level throughout the remainder of the test. 


92 


/ 


Test Conditions 


Tests were conducted under the following conditions. 


Condition 

Sink Speed 
m/sec (ft/sec) 

Fully-extended 
Gear Pressure 
kPa (lb/in^) 

Mode 

1. 

0. 305 (1) 

1048 (152) 

Passive 

2. 

0. 914 (3) 

1048 (152) 

Passive 

3. 

1. 524 (5) 

1048 (152) 

Passive 

4. 

0. 305 (1) 

1048 (152) 

Active 

5. 

0. 914 (3) 

1048 (152) 

Active 

6. 

1.524 (5) 

1048 (152) 

Active 

7. 

0. 305 (1) 

1931 (280) 

Passive 

8. 

0. 914 (3) 

1931 (280) 

Passive 

9. 

1.524 (5) 

1931 (280) 

Passive 

10. 

0.305 (1) 

1931 (280) 

Active 

) 

11. 

0. 914 (3) 

1931 (280) 

Active 

12. 

1.524 (5) 

1931 (280) 

Active 


Test Procedure 

The procedure for conducting the drop tests is given in Appendix C. 

Test Results 

The peak forces produced and percentage reduction due to active control 
are presented in Table VL Comparative plots of force, strut deflection, 
and pressure are shown in Figures 47 through 64. The oscillograph 
recordings from which those plots were constructed are included as 
Appendix D . 
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TABLE VL - PEAK FORCES AND PERCENTAGE REDUCTION 
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Figure 47. Net Force, Conditions 1, 
(From W/G Accel #1) 













Figure 48. Net Force, Conditions 2, 
(From W/G Accel #1) 
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Figure 49. Net Force, Conditions 
(From W/G Accel #1) 
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Figure 50„ Net Force, Conditions 
(From W/G Accel #1) 
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Figure 51. Net Force, Conditions 8, 11 
(From W/G Accel #1) 
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Figure 52. Net Force, Conditions 9, 12 
(From W/G Accel #1) 
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Figure 56. Gear Hydraulic Pressure, Conditions 
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Figure 59. Strut Deflection, Conditions 1 
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Figure 60. Strut Deflection, Conditions 2 
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Figure 62. Strut Deflection, Conditions 
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Figure 63. Strut Deflection, Conditions 8, 11 




Figure 64. Strut Deflection, Conditions 9, 12 


The apparent lack of reduction at the second peak under conditions 3 and 6 
is due to the fact that a resonance at approximately 20 Hz is evident at 
these conditions, which masks the performance of the gear both passively 
and when actively controlled. Reference 2 indicates that a resonance of 
20 Hz is apparent in the fore and aft direction when the gear is dropped 
on an inclined surface. This resonance can be suppressed by inclusion 
of appropriate compensation networks in the controller. 
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CONCLUSIONS 


This report has presented the analysis of an active control 
landing gear, the detailed design of an electronic controller to produce the 
active control, and the results of tests of the complete system. The results 
indicate that the force reduction varies from 9 to 31% depending on the air- 
craft sink speed and the static gear pressure. 

It is apparent that the effectiveness of the controller is highly dependent on 
the dy nam ic compensation employed; and, such compensation must be 
optimized for the particular aircraft landing gear system in which the 
controller is used. The compensation must be such that the dynamic re- 
sponse of the controller is quite high to control the impact force while still 

providing stable operation. 
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Appendix A 

MICROPROCESSOR PROGRAM LISTING 


ASSN 0000 8000 
CU 0082 

EB 00EP 

PORTS OOE? 

P0RT6 00EA 
OUT I 0040 
0UT2 OOA1 
0UT3 OOP 0 
OUT4 OOP1 
OUTS OOP2 
OUT6 0093 
IOBAS OOAO 
NPASE 8000 
N10U 0000 
NHIGH 0080 
NSTAT OOA1 


NFLA6 

OOA7 

NBIAS 

0O7F 

BASE 

F700 

SCA 

F700 

FCR 

F701 

LCHAN 

F702 

CLR 

F703 

ADDAT 

F704 

DACO 

F708 

BAC1 

F70 A 

NUXO 

0001 

NUX1 

0000 

NUX2 

0002 

NUX3 

0003 

XHAX 

OPBP 

XTHR 

0038 

TVO 

0085 

TV1 

OOC8 

TV2 

002F 


TVEXP 0036 
STACK 3FFF 
RAN 3F80 
ACCEL 8000 
FLIN 3F80 
U6VEL 3F84 
STRUT 3F86 
SINK 3F88 
BXTHR 3F8A 
BXNAX 3F8C 


Assembled at 0000 


PE 

3F8E 

TRANS 

3F92 

START 

OOOO 

LI 

OOIE 

L8 

0O4B 

L9 

0088 

L10 

0096 

L1 1 

0099 

L4 

OODE 

L5 

0116 

L6 

0123 

L! 

0127 

L7 

0131 

IN3 

0143 

Ml 

0 1 4 C 

M2 

0165 

INI 

0174 

M3 

017B 

SUB2 

0189 

SHL 

0190 

MATH 

0195 

UAIT 

0197 

P6A 

01 AC 

R2 

01B6 


0000 

0000 

0000 

0000 

oooo 

0000 

0000 

0000 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 

oooo 


0010 

t 









0020 

* DEFINE CONSTANTS 







0030 

* 









0040 

* 


I/O PORTS 







0050 

CU 

EQU 

82H 

; GROUP 

2 CONTROL 

WORD 



0060 

EB 

EQU 

OEBH 

; GROUP 

2 CONTROL 

UORD 

ADDR. 


0070 

P0RT5 

EQU 

0E9H 

; PORT 

5 ADDR 





0080 

P0RT6 

EQU 

OEAH 

; PORT 

6 ADDR 





0090 

* 









0100 

* 

PORT 

6 OUTPUTS 

CTLR 

CTLR 

T/D 

LAND 


01 to 

* 



ENA 

RESET 

LAMP LAMP INT.ENA 

SL .ENA 

0120 

OUT 1 

EQU 

40H 

; OFF 

ON 

OFF 

OFF 

DIS 

DIS 

0130 

0UT2 

EQU 

0A1H 

; ON 

OFF 

ON 

OFF 

DIS 

ENA 

0140 

0UT3 

EQU 

90H 

; ON 

OFF 

OFF 

ON 

DIS 

DIS 

0150 

OUT 4 

EQU 

91 H 

; ON 

OFF 

OFF 

ON 

DIS 

ENA 

0160 

0UT5 

EQU 

92H 

; ON 

OFF 

OFF 

ON 

ENA 

DIS 

0170 

0UT6 

EQU 

93H 

; ON 

OFF 

OFF 

ON 

ENA 

ENA 

0180 

* 









0190 

* 

HATH 

BOARD- SC 

310 PARAMETERS 





0200 

IOBAS 

EQU 

OAOH 

; I/O BASE ADDR 




0210 

XBASE 

EQU 

8000H 

; MEMORY BASE 

ADDR 



0220 

HLOU 

EQU 

OOH 

; LS BYTE MEMORY 

BASE 

ADDR 


0230 

MHIGH 

EQU 

80H 

; HS BYTE MEMORY 

BASE 

ADDR 


0240 

HSTAT 

EQU 

0A1H 

; STATUS BYTE 





0250 

HFLAG 

EQU 

0A7H 

: FLAG 

BYTE 





0260 

HBIAS 

EQU 

7FH 

; FLOATING POINT 

EXPONENT BIAS 


0270 

* 










116 


oooo 


0280 * 

ANALOG I/O BOARD 

- 

SBC 732 PARAMETERS 

0000 


0290 BASE 

EQU 

0F700H 


MEMORY BASE ADDRESS 

0000 


0300 SCA 

EQU 

BASE+O 


A/D COMMAND-STATUS REGISTER 

0000 


0310 FOR 

EQU 

BASE + 1 


MUX ADDR AND GAIN REGISTER 

oooo 


0320 L CHAN 

EQU 

BASE + 2 


LAST CHANNEL REGISTER 

oooo 


0330 CLR 

EQU 

BASE+3 


CLEAR INTERRUPTS 

oooo 


0340 ADDAT 

EQU 

BASE + 4 


A/D DATA REGISTER 

oooo 


0350 DACO 

EQU 

BASE +8 


DACO OUTPUT 

oooo 


0360 DAC1 

EQU 

BASE+10 


DAC1 OUTPUT 

oooo 


0370 HUXO 

EQU 

01 


A/D GAIN=1 , POINT MUX TO U/G ACCEL 

oooo 


0380 MUX 1 

EQU 

00 


A/D GAIN=1 , POINT MUX TO U/G VEL 

oooo 


0390 MUX2 

EQU 

02 


A/D GAIN=1 , POINT MUX TO STRUT POS 

oooo 


0400 MUX3 

EQU 

03 


A/D GAIN=1, POINT MUX TO SINK RATE 

oooo 


0410 * 





oooo 


0420 * 

STRUT PARAMETERS 



oooo 


0430 XHAX 

EQU 

09B9H 

m 

f 

MAX STROKE= 8.5 IN. = 6.078 VOLTS 

oooo 


0440 XTHR 

EQU 

003BH 

m 

f 

THRESHOLD= 0.2 IN. = 0.143 VOLTS 

oooo 


0450 * 





oooo 


0460 * 

TRANSITION VELOCITY SCALE FACTOR 

oooo 


0470 * 

=0.00004191 DECIMAL 

oooo 


0480 * 

= 1.010111 1 1100100010000101 * 2=* : *-1 5 BINARY 

oooo 


0490 TOO 

EQU 

S5H 

f 

LS BYTE (H) 

oooo 


0500 TV1 

EQU 

0C8H 

9 

9 

LS BYTEM (M+1 ) 

oooo 


0510 TV2 

EQU 

2FH 

9 

MS BYTE (M+2) 

oooo 


0520 TVEXP 

EQU 

36H 

9 

EXPONENT (M+3) 

oooo 


0530 * 





oooo 


0540 * 

RAM 

MEMORY ASSIGNMENTS 

oooo 


0550 STACK 

EQU 

3FFFH 

■ 

9 

INITIAL STACK POINTER 

oooo 


0560 RAM 

EQU 

3F80H 

9 

START OF RAM SCRATCH AREA 

oooo 


0570 ACCEL 

EQU 

MBASE 

• 

r 

U/G ACCEL 

oooo 


0580 FLIM 

EQU 

RAM 

f 

LIMIT FORCE COMMAND 

oooo 


0590 WGVEL 

EQU 

RAM+4 

w 

f 

U/G VELOCITY 

oooo 


0600 STRUT 

EQU 

RAH+6 

f 

STRUT POSITION 

oooo 


0610 SINK 

EQU 

RAM + 8 

9 

INITIAL SINK RATE 

oooo 


0620 BXTHR 

EQU 

RAH+10 

• 

9 

STRUT THRESHOLD * 16 

oooo 


0630 BXMAX 

EQU 

RAM+12 

■ 

9 

MAX STRUT STROKE * 16 

oooo 


0640 PE 

EQU 

RAM+14 

9 

POTENTIAL ENERGY- 4 BYTES 

oooo 


0650 TRANS 

EQU 

RAH+18 

■ 

r 

TRANSITION VELOCITY- 4 BYTES 

oooo 


0660 * 





oooo 


0670 * START OF 

MAIN PROGRAM 



oooo 


0680 * 





OOOO F3 


0690 START 

DI 

.'DISABLE INTERRUPTS 

0001 21 

FF 3F 

0700 

LXI 

H, STACK : INIT 

STACK POINTER 

0004 F9 


0710 

SPHL 




0005 


0720 * 


CONFIGURE I/O PORTS 

0005 3E 

82 

0730 

MVI 

A , CU 



0007 D3 

EB 

0740 

OUT 

EB 



0009 


0750 * 


INIT 

MATH BOARD 

0009 3E 

00 

0760 

MVI 

A,MLOU ; SET 

MEMORY BASE ADDR 

OOOB D3 

A1 

0770 

OUT 

MSTAT 



OOOO 3E 

80 

0780 

MVI 

A,MHIGH 





OOOF H3 A2 


07?0 

OUT 

HSTAT+1 



0011 21 38 

00 

0800 

LXI 

H,XTHR 

; HULT STRUT THRESHOLD BY 16 


0014 CD 90 

01 

0810 

CALL 

SHL 

; FOR LATER USE 


0017 22 BA 

3F 

0820 

SHLD 

BXTHR 



001 A 3E 40 


0830 

HVI 

A, 0UT1 

; SET LIGHTS, SUITCHES 


OOIC D3 EA 


0840 

OUT 

P0RT6 



OOIE DB E? 


0850 LI 

IN 

P0RT5 

{CONTROLLER ENABLED? 


0020 IF 


0860 

RAR 




0021 D2 IE 

00 

0870 

JNC 

LI 

; NO , KEEP LOOKING 


0024 


0880 * 





0024 


0890 * 



CONTROLLER HAS BEEN ENABLED 


0024 


0?00 * 





0024 3E 02 


0?10 

HVI 

A,MUX2 

; YES, GET STRUT POSITION FOR 


0026 CD 74 

Ot 

0920 

CALL 

INI 

{LANDINfl/TAKEOFF DETERMINATION 


002? 2A 8A 

3F 

0930 

LHLD 

BXTHR 

{GET STRUT THRESHOLD 


O02C ED 


0940 

XCHG 


{PUT IN DE 


O02D 2A 86 

3F 

0950 

LHLD 

STRUT 

{LOAD HL UITH STRUT POSN 


0030 CD 8? 

01 

0960 

CALL 

SUB2 

{CALC: THRESHOLD - STRUT 


0033 DA 27 

01 

0970 

JC 

L2 

{TAKING OFF 


0036 


0980 * 





0036 


0990 * 



LANDING, HAKE PREPARATIONS 


0036 


1000 * 





0036 3E 03 


1010 

HVI 

A, HUX3 

{GET INITIAL SINK RATE 


0038 CD 74 

01 

1020 

CALL 

INI 



O03B 22 88 

3F 

1030 

SHLD 

SINK 

{STORE IT 


O03E 21 B? 

0? 

1040 

LXI 

H,XMAX 

{HULT XHAX BY 16 TO SHIFT INTO UPPER 12 BITS 

0041 CD ?0 

01 

1050 

CALL 

SHL 



0044 22 8C 

3F 

1060 

SHLD 

BXHAX 

{STORE IT 


0047 


1070 * 





0047 


1080 * 



ENABLE INTEGRATOR, START 


0047 


1090 * 



ENERGY CALCULATIONS. 


0047 


1100 * 





0047 3E ?2 


1110 

HVI 

A , 0UT5 

{ENABLE INT 


004? D3 EA 


1120 

OUT 

P0RT6 



O04B CD 43 

01 

1130 L8 

CALL 

IN3 

{GET ACCEL, 11/ G VEL, STRUT POSN 

FROM A/D 

O04E EB 


1140 

XCHG 


{CALC POTENTIAL ENERGY. SAVE HL 

IN DE 

004F 2A 8C 

3F 

1150 

LHLD 

BXMAX 

{GET HAX STROKE 


0052 EB 


1160 

XCHG 


. {PUT IN DE 


0053 7B 


1170 

HOV 

A,E 

{CALC: XMAX - STRUT POSN 


0054 ?5 


1180 

SUB 

L 



0055 6F 


1190 

HOV 

M 



0056 7A 


1200 

HOV 

A,D 



0057 ?C 


1210 

SBB 

H 



0058 67 


1220 

HOV 

H, A 



005? 22 04 

80 

1230 

SHLD 

MBASE+4 ; STORE IN HATH BOARD 


O05C 


1240 *. 



ACCEL IS ALREADY IN MATH BOARD AT HBASE+0,1 

O05C AF 


1250 

XRA 

A 

{MULTIPLY 


O05D CD ?5 

01 

1260 

CALL 

HATH 

{NOU HAVE PE AS A 32-BIT UORB 


0060 2A 00 

BO 

1270 

LHLD 

HBASE 

{SAVE IT IN RAH 


0063 22 8E 

3F 

1280 

SHLD 

PE 



0066 2A 02 

80 

1290 

LHLD 

HBASE+ 

2 
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JCALC KINETIC ENERGY 


0069 22 90 3F 
O06C 2A 88 3F 
O06F EB 

0070 2A 84 3F 
0073 CD 89 01 
0076 22 00 80 

0079 22 04 80 
007 C AF 

O07D CD 95 01 

0080 
0080 
0080 

0080 06 03 
0082 21 91 3F 
0085 11 03 80 

0088 1 A 

0089 BE 

O08A C2 96 00 
O08D IB 
008E 2B 
O08F 05 

0090 C2 88 00 
0093 C3 99 00 
0096 D2 4B 00 
0099 

0099 

0099 

0099 2A 80 3 F 

009C 22 08 F7 

009F 3E 93 

OOA1 D3 EA 

OOA3 

OOA3 

OOA3 

OOA3 

OOA3 2A 80 3F 
OOA6 22 00 80 
OOA9 21 00 00 
OOAC 22 02 80 
OOAF 3E 08 
OOB1 CD 95 01 
OOB4 3E 06 
OOB6 CD 95 01 
OOB9 3E 85 
OOBB 32 04 80 
OOBE 3E C8 
OOCO 32 05 80 
OOC3 3E 2F 
OOC5 32 06 80 
OOCB 3E 36 
OOCA 32 07 80 


1300 

SHLD 

PE+2 

1310 

LHLD 

SINK 

1320 

XCHG 


1330 

LHLD 

UGVEL 

1340 

CALL 

SUB2 

1350 

SHLD 

MBASE 

1360 

SHLD 

MBASE+4 

1370 

XRA 

A 

1380 

CALL 

MATH 

1390 * 



1400 * 



1410 * 



1420 

MVI 

B , 3 

1430 

LXI 

H, PE+3 

1440 

LXI D 

, HBASE+3 

1450 L9 

LDAX 

D 

1460 

CMP 

M 

1470 

JNZ 

LI 0 

1480 

DCX 

D 

1490 

DCX 

H 

1500 

DCR 

B 

1510 

JNZ 

L9 

1520 

JMP 

L 1 1 

1530 LI 0 

JNC 

L8 

1540 * 



1550 * 



1560 * 



1570 L1 1 

LHLD 

FLIM 

1580 

SHLD 

DACO 

1590 

MVI 

A,0UT6 

1600 

OUT 

P0RT6 


{CALC: SINK RATE - U/G VEL 
r LOAD INTO MATH BOARD, TWO PLACES 

{MULT 

:nou have ke at mbase+o, 1 , 2,3 

DO A BYTE-BY-BYTE COMPARISON TO TEST IF 
PE GREATER THAN KE. DO NOT BOTHER TO TEST 
LS BYTE, IT HAS NO USEFUL DATA. 

{SET A BYTE COUNTER 
, ’POINT HL TO PE MSB 
JPOINT DE TO KE MSB 
,-PE GREATER THAN KE? 


;trv another byte 

{ TESTED 3 BYTES? 

{NO, LOOP BACK 

;G0 GET NEW INPUTS, TRY AGAIN 
TIME TO INITIATE ACTIVE CONTROL 
; ISSUE LIMIT FORCE COMMAND 
, -ENABLE SERVOLOOP 


1610 * 


1620 * GEAR IS NOU UNDER ACTIVE CONTROL 
1630 * 


1640 * 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 


CALC TRANSITION VELOCITY 
LHLD FLIM {GET LIMIT FORCE CMD ( FL I > 

SHLD MBASE ; LOAD INTO MATH BOARD 

LXI H.O 
SHLD MBASE+2 

NVI A, 8 {CONVERT TO FLOATING POINT 
CALL MATH 

MVI A, 6 ; SQUARE IT 

CALL MATH 

MI A, TVO {LOAD TRANSITION VELOCITY SCALE 

STA MBASE+4 {FACTOR INTO NATH BOARD 

MVI A, TV1 

STA MBASE+5 

MVI A , TV2 

STA MBASE+6 

MVI A,TVEXP 

STA MBASE+7 
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OQCB 3E 02 
OOCF CD ?5 01 
O0D2 2A 00 80 
O0D5 22 92 3F 
OOD8 2A 02 80 
OODB 22 94 
OODE 
OODE 
OODE 

OODE 2A 88 3F 
OOE1 EB 
OOE2 3E 00 
OOE4 CD 74 01 
OOE7 CD 89 01 
OOEA 22 00 80 
OOED 21 00 00 
OOFO 22 02 80 
OOF3 3E 08 
O0F5 CD 95 01 
OOF8 2A 92 3F 
OOFB 22 04 80 
OOFE 2A 94 3F 
0101 22 06 80 
0104 3E OA 
0106 CD 95 01 
0109 DB A1 
010B £6 20 
010D CA DE 00 



1810 

HVI 

A, 2 ; MULT ! 

01 

1820 

CALL 

HATH 

80 

1830 

LHLD 

HBASE ; STORE 

3F 

1840 

SHLD 

TRANS 

80 

1850 

LHLD 

MBASE+2 

3F 

1860 

SHLD 

TRANS+2 


1870 

* NOU HAVE 

TRANSITION VEL 


1880 * 
1890 * 
1900 14 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 


bu D 1 H*\ * lunrmunw • " * ~ 

FOR DETERMINING START OF TRANSITION. 


U/G VEL) 


LHLD 

XCHG 

HVI 

CALL 

CALL 

SHLD 

LXI 

SHLD 

HVI 

CALL 

LHLD 

SHLD 

LHLD 

SHLD 

HVI 

CALL 

IN 

ANI 

JZ 


SINK 

A , MUX 1 

INI 

SUB2 

HBASE 

H,0 

HBASE+2 

A, 8 

MATH 

TRANS 

HBASE+4 

TRANS+2 

MBASE+6 

A,OAH 

HATH 

MSTAT 

20H 

L4 


0110 


2090 * 



0110 


2100 * TRANSITION 

PHASE 

0110 


2110 * 



0110 2A 

80 3F 

2120 

LHLD 

FLIH 

0113 11 

FD FF 

2130 

LXI 

D,-3 

0116 22 

08 F7 

2140 L5 

SHLD 

DACO 

0119 19 


2150 

DAD 

D 

Oil A DA 

16 01 

2160 

JC 

L5 

01 ID 21 

00 00 

2170 

LXI 

H,0 

0120 22 

08 F7 

2180 

SHLD 

DACO 

0123 00 


2190 L6 

NOP 


0124 C3 

23 01 

2200 

JHP 

L6 

0127 


2210 * 



0127 


2220 * TAKEOFF HODE 

0127 


2230 * 



0127 21 

00 00 

2240 L2 

LXI 

H ,0 

0 1 2 A 22 

08 F7 

2250 

SHLD 

DACO 

01 2D 3E 

A1 

2260 

HVI 

A,0UT2 

012F D3 

EA 

2270 

OUT 

P0RT6 

0131 3E 

02 

2280 L7 

HVI 

A, HUX2 

0133 CD 

74 01 

2290 

CALL 

INI 

0136 EB 


2300 

XCHG 


0137 2A 

8A 3F 

2310 

LHLD 

BXTHR 


; GET SINK RATE 
; SAVE IN DE 
; GET U/G VELOCITY 

;CALC: SINK RATE - U/G VEL 
; CONVERT TO FLOATING POINT 


;ITS NOU IN MBASE+O, 1 ,2,3 
; LOAD TRANS VEL INTO HATH BOARD 


; COMPARE AGAINST (SINK-U/G VEL) 

:GET STATUS, MASK 'LESS THAN' BIT 
; IS (SINK-U/G VEL) .LT. TRANS VEL? 
:N0, CONTINUE LOOPING 
YES, TIME TO START TRANSITION 


; LOAD LIMIT FORCE COMMAND 

; SET RAHP RATE 

; OUTPUT CMD TO DAC 

DECREASE LIHIT FORCE CMD 

; LOOP UNTIL CMD = 0 

jSET LIHIT FORCE CMD EXACTLY = 0 

: STAY I.N A LOOP UNTIL A RESET OCCURS 


; COWHAND A ZERO LIHIT FORCE 

; ENABLE SERVO LOOP, LEAVE ENABLED 
; UNTIL STRUT POSITION LESS THAN THRESHOLD 
•GET STRUT POSITION 

;PUT IN DE 

; LOAD HL UITH THRESHOLD 
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013A CD 8? 01 

013D D2 31 01 

0140 C3 00 00 

0143 

0143 

0143 

0143 

0143 3E 01 
0145 21 01 F7 
0148 77 
014? 2B 
0 1 4 A 36 01 
014C 7E 
01 4D 07 
014E D2 4C 01 
0151 36 00 
0153 2A 04 F7 
0156 22 00 80 
015? 22 80 3F 
015C 3E 00 
OISE 21 01 F7 

0161 77 

0162 2B 

0163 36 01 

0165 7E 

0166 07 

0167 D2 65 01 
016A 36 00 
016C 2A 04 F7 
0I6F 22 84 3F 
0172 3E 02 
0174 21 01 F7 

0177 77 

0178 2B 
017? 36 01 
0I7B 7E 
0I7C 07 

017D D2 7B 01 
0180 36 00 
0182 2A 04 F7 
0185 22 86 3F 
0188 C? 

018? 

018? 

018? 

018? 7B 
018A ?5 
018B 6F 
018C 7A 
018D ?C 
018E 67 


2320 

CALL 

SUB2 

2330 

JNC 

L7 

2340 
2350 * 

JMP 

START 


J CALC s STRUT - THRESHOLD 

{LOOP UNTIL STRUT EXTENDED FULLY 

JHAVE LIFTOFF, TURN OFF CONTROLLER 


2360 * ROUTINE TO INPUT AND STORE DATA FROM 
2370 * THREE HUX CHANNELS 
2380 * 


2 3 ?0 IN3 

MV I 

A, MUXO 

POINT MUX TO U/G ACCEL 

2400 

LXI 

H,FCR 

;POINT HL TO MUX/GAIN REGISTER 

2410 

MOV 

M, A 

;load register 

2420 

DCX 

H 

JPOINT HL TO CMD/STATUS REGISTER 

2430 

MVI 

M,01 

J START CONVERSION 

2440 Ml 

MOV 

A,M 

{READ STATUS 

2450 

RLC 


; DONE? 

2460 

JNC 

Ml 

;no, KEEP LOOPING 

2470 

MVI 

M,0 

; YES, RESET CONVERSION ENABLE 

2480 

LHLD 

ADDAT 

;get data 

24 ?0 

SHLD 

ACCEL 

;STO U/G ACCEL IN MATH BOARD, 

2500 

SHLD 

FLIM 

{ ALSO IN RAH 

2510 

MVI 

A , MUX 1 

{REPEAT FOR U/G VELOCITY 

2520 

LXI 

H,FCR 


2530 

MOV 

M, A 


2540 

DCX 

H 


2550 

MVI 

M,01 


2560 M2 

MOV 

A,M 


2570 

RLC 



2580 

JNC 

M2 


25?0 

MVI 

M, 0 


2600 

LHLD 

ADDAT 


2610 

SHLD 

UGVEL 

{STORE U/G VEL 

2620 

MVI 

A , MUX2 

{REPEAT FOR STRUT POSITION 

2630 INI 

LXI 

H,FCR 


2640 

MOV 

M, A 


2650 

DCX 

H 


2660 

MVI 

M, 01 


2670 M3 

MOV 

A ,H 


2680 

RLC 



26?0 

JNC 

M3 


2700. 

MVI 

M, 0 


2710 

LHLD 

ADDAT 


2720 

SHLD 

STRUT 


2730 

RET 



2740 * 




2750 * DOUBLE PRECISION 

SUBTRACT ROUTINE 

2760 * HL 

=DE-HL 



2770 SU62 

MOV 

A,E 


2780 

SUB 

L 


27?0 

MOV 

L r A 


2800 

MOV 

A,D 


2810 

SBB 

H 


2820 

MOV 

H, A 




01 8F C9 
0190 
0190 
0190 

0190 29 

0191 29 

0192 29 

0193 29 

0194 C9 

0195 
0195 
0195 
0195 
0195 03 
0197 DB 
0199 E6 
019B C2 
019E C9 
019F 
019F 
019F 
019F 
019F 
019F 
019F F3 
01 AO 3E 
01A2 CD 
01A5 CF 
01 A6 00 


AO 

A7 

01 

97 01 


00 

74 01 


2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 


RET 


* ROUTINE TO SHIFT VALUE IN HL LEFT 4 PLACES. 

* 

SHL DAD H 
H 
H 
H 


DAD 

DAD 

DAD 

DAD 

RET 


* ROUTINE TO ACTIVATE MATH BOARD AND UAIT FOR RESULT. 


2940 

* ACCUH HAS 

OPCODE. 


2950 

* 



; COMMAND MATH BOARD TO START 

2960 

HATH 

OUT 

IOBAS 

2970 

UAIT 

IN 

I OB AS +7 

; GET FLAG BYTE 

2980 


ANI 

01 

; CHECK BUSY BIT 

2990 


JNZ 

UAIT 

JSTAY IN LOOP UNTIL NOT BUSY 

3000 


RET 




3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 


* 

* 

* SPECIAL CHECK-OUT ROUTINES 

* 

* ROUTINE TO INPUT A VALUE FROM A/D, STORE IN RAN. 


DI 
MV I 
CALL 
RST 
NOP 


A, 00 
INI 
1 


; SELECT CHAN 0 


01A7 

00 



3120 


NOP 

DO PGA TEST ON A/D 

01A8 




3130 

* ROUTINE TO 

01A8 

F3 



3140 


DI 


01A9 

21 

01 

F7 

3150 


LXI 

H , FCR 

01 AC 

36 

00 


3160 

PGA 

MVI 

M,O0 

01 AE 

36 

CO 


3170 


HVI 

H,OCOH 

01B0 

C3 

AC 

01 

3180 


JHP 

PGA 

0IB3 

00 



3190 


NOP 

TO OUTPUT A VALUE TO DACO, DAC1 

01B4 



3200 

* 

ROUTINE 

0IB4 

F3 



3210 


DI 


01B5 

00 



3220 


NOP 


01B6 

21 

00 

00 

3230 

R2 

LXI 

H,0 

01B9 

22 

08 

F7 

3240 


SHLU 

DACO 

01BC 

22 

OA 

F7 

3250 


SHLD 

DAC1 

01BF 

00 



3260 


NOP 


OICO 

00 



3270 


NOP 


01C1 

00 



3280 


NOP 


01C2 

C3 

B6 

01 

3290 


JHP 

R? 
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Assembled at 3D 10 


ASSN 

3D1 0 8D10 

CU 

0082 

EB 

OOEB 

P0RT5 

OOE9 

P0RT6 

OOEA 

OUT! 

0040 

0UT2 

OOA1 

0UT3 

0090 

0UT4 

0091 

OUTS 

0092 

0UT6 

0093 

I OB AS 

OOAO 

NBASE 

8000 

NLOU 

OOOO 

MHIGH 

0080 

NSTAT 

OOA1 

MFLAG 

OOA7 

MB I AS 

007F 

BASE 

F70O 

SCA 

F70O 

FCR 

F701 

LCHAN 

F702 

CLR 

F703 

ABDAT 

F704 

DACO 

F708 

DACt 

F70A 

MUXO 

000 1 

MUX t 

0000 

NUX2 

0002 

MUX3 

0003 

XHAX 

09B9 

XTHR 

0038 

TVO 

0085 

TV1 

OOC8 

TV2 

002F 

TVEXP 

0036 

STACK 

3FFF 

RAH 

3F80 

ACCEL 

8000 

FLIM 

3F80 

U6VEL 

3F84 

STRUT 3F86 

SINK 

3F88 

BXTHR 

3F8A 

BXHAX 

3F8C 

PE 

3F8E 

TRANS 

3F92 

START 

3D10 

L! 

3D2E 

LB 

3D5B 

L9 

3B98 
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LI 0 

3DA6 

LI1 

3DA? 

L4 

3DEE 

LS 

3E26 

L4 

3E33 

L2 

3E37 

L7 

3E41 

IN3 

3E53 

Ml 

3ESC 

M2 

3E75 

INI 

3E84 

M3 

3E8B 

SUB2 

3E99 

SHL 

3EA0 

MATH 

3EA5 

UAIT- 

3EA7 

PGA 

3EBC 

R2 

3EC6 

3910 


3910 


3910 



3D10 

3D10 

3D1 0 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 

3910 


0010 * 

0020 * DEFINE CONSTANTS 


0030 

* 



0040 

4 


I/O PORTS 

0050 

cu 

EQU 

82H 

0060 

EB 

EQU 

OEBH 

0070 

PORTS 

EQU 

0E9H 

0080 

P0RT6 

EQU 

OEAH 

0070 

* 



0100 

* 

PORT 

6 OUTPUTS 

0110 

* 



0120 

0UT1 

EQU 

40H 

0130 

0UT2 

EQU 

0A1H 

0140 

0UT3 

EQU 

90H 

0150 

0UT4 

EQU 

91 H 

0160 

0UT5 

EQU 

92H 

0170 

0UT6 

EQU 

93H 


GROUP 2 CONTROL UORH 
GROUP 2 CONTROL WORD ADDR. 
PORT 5 ADDR 
PORT 6 ADDR 


CTLR 

CTLR 

T/0 

ENA 

RESET 

LAMP 

; OFF 

ON 

OFF 

; ON 

OFF 

ON 

; ON 

OFF 

OFF 

; ON 

OFF 

OFF 

; ON 

OFF 

OFF 

: ON 

OFF 

OFF 


LAND 

LAMP 

INT.ENA 

SL .ENA 

OFF 

DIS 

DIS 

OFF 

DIS 

ENA 

ON 

DIS 

DIS 

ON 

DIS 

ENA 

ON 

ENA 

DIS 

ON 

ENA 

ENA 


0100 

* 



0190 

* 

MATH 

BOARD- SC 310 

0200 

IOBAS 

EQU 

OAOH 

0210 

MBASE 

EQU 

8000H 

0220 

MLOU 

EQU 

OOH 

0230 

MHIGH 

EQU 

80H 

0240 

MSTAT 

EQU 

0A1H 

0250 

MFLAG 

EQU 

0A7H 

0260 

MBIAS 

EQU 

7FH 

0270 

* 



0280 

* 

ANALOG I/O BOARD- 

0290 

BASE 

EQU 

0F700H 

0300 

SCA 

EQU 

BASE+O 

0310 

FCR 

EQU 

BASE+1 

0320 

LCHAN 

EQU 

BASE+2 

0330 

CLR 

EQU 

BASE+3 


PARAMETERS 

I/O BA5E ADDR 

HEHORY BASE ADDR 

LS BYTE MEMORY BASE ADDR 

MS BYTE MEMORY BASE ADDR 

STATUS BYTE 

FLAG BYTE 

FLOATING POINT EXPONENT BIAS 

SBC 732 PARAMETERS 
MEMORY BASE ADDRESS 
A/D COMMAND-STATUS REGISTER 
MUX ADDR AND GAIN REGISTER 
LAST CHANNEL REGISTER 
CLEAR INTERRUPTS 
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0340 

0350 

0360 

0370 

0380 

0390 

0400 

0410 

0420 

0430 

0440 

0450 

0460 

0470 

0480 

0490 


ADD A T 
DACO 
DAC1 
MUXO 
MUX1 
HUX2 
MUX3 
♦ 

* 

XMAX 
X T HR 
* 

* 

* 

* 

TOO 


EQU BASE+4 
EQU PASE+8 
EQU BASE + 1 0 
EQU 01 
EQU 00 
EQU 02 
EQU 03 

STRUT PARAMETERS 
EQU 09B9H 
EQU 0038H 


A/D DATA REGISTER 
DACO OUTPUT 
»AC1 OUTPUT 
A/D 6AIN=1 , POINT 
A/D GAIN= 1 , POINT 
A/D GAIN=1 , POINT 
A/D GA I N = 1 , POINT 


POINT MUX TO U/G ACCEL 
POINT NUX TO U/G VEL 
POINT MUX TO STRUT POS 
POINT MUX TO SINK RATE 


MAX STROKE= 8.5 IN. = 
THRESHOLB= 0.2 IN. = < 


6.078 VOLTS 
'.143 VOLTS 


3D1 0 



0500 

TV1 

EQU 

0C8H 

3D10 



0510 

TV2 

EQU 

2FH 

3010 



0520 

TVEXP 

EQU 

36H 

3010 



0530 

* 



3010 



0540 

* 

RAM 

MEMORY A 

3010 



0550 

STACK 

EQU 

3FFFH 

3010 



0560 

RAH 

EQU 

3F80H 

3010 



0570 

ACCEL 

EQU 

MBASE 

3010 



0580 

FLIM 

EQU 

RAM 

3010 



0590 

UGVEL 

EQU 

RAM+4 

3010 



0600 

STRUT 

EQU 

RAM+6 

3010 



0610 

SINK 

EQU 

RAM + 8 

3010 



0620 

BXTHR 

EQU 

RAM+10 

3010 



0630 

BXMAX 

EQU 

RAH+12 

3010 



0640 

PE 

EQU 

RAM+14 

3010 



0650 

TRANS 

EQU 

RAH+18 

3010 



0660 

* 



3010 



0670 

* START OF 

MAIN PROI 

3010 



0680 

* 



3010 F3 



0690 

START 

BI 


3011 21 

FF 

3F 

0700 


LX I 

H, STACK 

3014 F9 



0710 


SPHL 


3015 



0720 

* 



3015 3E 

82 


0730 


MVI 

A , CU 

3017 D3 

EB 


0740 


OUT 

ED 

3019 



0750 

* 



3019 3E 

00 


0760 


MVI 

A , MLOU 

301 B D3 

A1 


0770 


OUT 

MSTAT 

3D 1 D 3E 

80 


0780 


MVI 

A,HHIGH 

3D1F D3 

A2 


0790 


OUT 

MSTAT+I 

3021 21 

38 

00 

0800 


LXI 

H, XTHR 

3024 CD 

AO 

3E 

0810 


CALL 

SHL 

3027 22 

8A 

3F 

0820 


SHLD 

BXTHR 

302A 3E 

40 


0830 


MVI 

A, OUT! 

3D2C 03 

EA 


0840 


OUT 

P0RT6 


TRANSITION VELOCITY SCALE FACTOR 
=0.00004191 DECIMAL 

= 1.010111 11100100010000101 * 2* * — 1 5 BINARY 
EQU 85H ; |_S BYTE (M) 

EQU 0C8H ; LS BYTE+1 <M+1) 

EQU 2FH ; MS BYTE (H+2) 

EQU 36H : EXPONENT (H+3) 


INITIAL STACK POINTER 
START OF RAM SCRATCH AREA 
U/G ACCEL 

LIMIT FORCE COMMAND 
U/G VELOCITY 
STRUT POSITION 
INITIAL SINK RATE 
STRUT THRESHOLD * 16 
MAX STRUT STROKE * 16 
POTENTIAL ENERGY- 4 BYTES 
TRANSITION VELOCITY- 4 BYTES 


I DISABLE INTERRUPTS 
JINIT STACK POINTER 

CONFIGURE I/O PORTS 


IN IT HATH BOARD 
[SET MEHORY BASE 


;hult strut threshold by h 

;FOR LATER USE 
JSET LIGHTS, SUITCHES 
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3D2E 

66 

o> 

UJ 


0850 LI 

IN 

PORTS 

; CONTROLLER ENABLED? 

3630 

IF 



0860 

RAR 



3D31 

62 

2E 

3D 

0870 

JNC 

LI 

;N0, KEEP LOOKING 

3D34 




0880 * 




3534 




0890 * 


CONTROLLER HAS BEEN ENABLED 

3534 




0900 * 




3534 

3E 

02 


0910 

MV I 

A ,MUX2 

; YES, GET STRUT POSITION FOR 

3536 

CD 

84 

3E 

0920 

CALL 

INI 

JLANDING/TAKEOFF DETERMINATION 

353? 

2A 

8A 

3F 

0930 

LHLD 

BXTHR 

;GET STRUT THRESHOLD 

353C 

EB 



0940 

XCHG 


; PUT IN DE 

353D 

2A 

86 

3F 

0950 

LHLD 

STRUT 

; LOAD HL WITH STRUT POSN 

3540 

CD 

9? 

3E 

0960 

CALL 

SUB2 

;CALC: THRESHOLD - STRUT 

3543 

DA 

37 

3E 

0970 

JC 

L2 

;TAKING OFF 

3546 




0980 * 




3546 




0990 * 



LANDING, HAKE PREPARATIONS 

3546 




1000 * 




3546 

3E 

03 


1010 

M VI 

A,HUX3 

;GET INITIAL SINK RATE 

3648 

CD 

84 

3E 

1020 

CALL 

INI 


3546 

22 

88 

3F 

1030 

SHLD 

SINK 

; STORE IT 

364E 

21 

B? 

0? 

1040 

LXI 

H, XMAX 

;MULT XMAX BY 16 TO SHIFT INTO UPPER 12 BITS 

3651 

CD 

AO 

3E 

1050 

CALL 

SHL 


3554 

22 

8C 

3F 

1060 

SHLD 

BXMAX 

; STORE IT 

3557 




1070 * 




3557 




1080 * 



ENABLE INTEGRATOR, START 

3557 




1090 * 



ENERGY CALCULATIONS. 

3657 




1100 * 




3657 

3E 

92 


1110 

HOI 

A,0UT5 

JENABLE INT 

3559 

63 

EA 


1120 

OUT 

P0RT6 


3D5B 

CD 

53 

3E 

1130 L8 

CALL 

IN3 

; SET ACCEL, U/G VEL, STRUT POSN FROM A/D 

365E 

EB 



1140 

XCHB 


; CALC POTENTIAL ENERGY. SAVE HL IN DE 

365F 

2A 

8C 

3F 

1150 

LHLD 

BXMAX 

; GET MAX STROKE 

3562 

EB 



1160 

XCHG 


;PUT IN DE 

3663 

7B 



1170 

MOV 

A,E 

; CALC : XHAX - STRUT POSN 

3664 

95 



1180 

SUB 

L 


3665 

6F 



1190 

MOV 

L,A 


3666 

7A 



1200 

MOV 

A,D 


3667 

?C 



1210 

SBB 

H 


3668 

67 



1220 

MOV 

H,A 


3669 

22 

04 

80 

1230 

SHLD 

MBASE+4 

JSTORE IN HATH BOARD 

366C 




1 240 * 



ACCEL IS ALREADY IN MATH BOARD AT HBASE+0,1 

356C 

AF 



1250 

XRA 

A 

^MULTIPLY 

3666 

CD 

A5 

3E 

1260 

CALL 

MATH 

;NOU HAVE PE AS A 32-BIT UORD 

3670 

2A 

00 

80 

1270 

LHLD 

MBASE 

; SAVE IT IN RAM 

3673 

22 

BE 

3F 

1280 

SHLD 

PE 


3576 

2A 

02 

80 

1290 

LHLD 

MBASE+2 


367? 

22 

90 

3F 

1300 

SHLD 

PE+2 


3570 

2A 

88 

3F 

1310 

LHLD 

SINK 

;CALC KINETIC ENERGY 

357F 

EB 



1320 

XCHG 



3580 

2A 

84 

3F 

1330 

LHLD 

UGVEL 


3683 

CD 

9? 

3E 

1340 

CALL 

SUB2 

;CALC: SINK RATE - U/G VEL 

3686 

22 

00 80 

1350 

SHLD 

MBASE 

; LOAD INTO MATH BOARD, TUO PLACES 
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308? 22 04 80 

3DBC AF 

3I8D CO A3 3E 

3090 

3090 

3090 

3D90 06 03 

3092 21 91 3F 

3093 11 03 80 

3098 1A 

3099 BE 

309A C2 A6 30 
3090 IB 
3B9E 2B 
309F 05 
30AO C2 98 30 
3DA3 C3 A9 30 
30A6 02 5B 30 
30A9 


1360 

SHLD 

MBASE+4 

1370 

XRA 

A 

1380 

CALL 

MATH 

1390 * 



1400 * 



1410 * 



1420 

MVI 

B,3 

143D 

LXI 

H,PE+3 

1440 

LXI 0 

,HBASE+3 

1450 L9 

LDAX 

0 

1460 

CMP 

M 

1470 

JNZ 

LI 0 

1480 

OCX 

D 

1490 

OCX 

H 

1300 

OCR 

B 

1510 

JNZ 

L9 

1520 

JMP 

L1 1 

1530 LtO 

JNC 

L8 


;mult 

;nou have ke at mbase+o, 1 , 2,3 

DO A BYTE-BY-BYTE COMPARISON TO TEST IF 
PE GREATER THAN KE. DO NOT BOTHER TO TEST 
LS BYTE, IT HAS NO USEFUL DATA. 

;SET A BYTE COUNTER 
;POINT HL TO PE MSB 
;POINT DE TO KE MSB 
JPE GREATER THAN KE? 


,-TRY ANOTHER BYTE 


1540 * 


30A9 

1550 * 


30A9 

1560 * 


30A9 2A 80 3F 

1370 Lit 

LHLD 

3DAC 22 08 F7 

1580 

SHLD 

30AF 3E 93 

1590 

MVI 

3DB1 03 EA 

1600 

OUT 

30B3 

161 0 * 


3DB3 

1620 * GEAR IS 1 

3DB3 

1630 * 


30B3 

1640 * 


30B3 2A BO 3F 

1650 

LHLD 

30B6 22 00 80 

1660 

SHLD 

30B9 21 00 00 

1670 

LXI 

3DBC 22 02 80 

1680 

SHLD 

3DBF 3E 08 

1690 

MVI 

30C1 CD A5 3E 

1700 

CALL 

3DC4 3E 06 

1710 

MVI 

30C6 CD A5 3E 

1720 

CALL 

30C9 3E 85 

1730 

MVI 

3DCB 32 04 80 

1740 

STA 

30CE 3E C8 

1750 

MVI 

3DD0 32 05 80 

1760 

STA 

30D3 3E 2F 

1770 

MVI 

3DD5 32 06 80 

1780 

STA 

30D8 3E 36 

1790 

MVI 

30DA 32 07 80 

1800 

STA 

30DD 3E 02 

1810 

MVI 

3DDF CD A5 3E 

1820 

CALL 

30E2 2A 00 80 

1830 

LHLD 

3DE5 22 92 3F 

1840 

SHLD 

30E8 2A 02 80 

1850 

LHLD 

30EB 22 94 3F 

1860 

SHLD 


FLIN 

DACO 

A,0UT6 

P0RT6 


; TESTED 3 BYTES? 

;no, loop back 

,-GO GET NEU INPUTS, TRY AGAIN 
TINE TO INITIATE ACTIVE CONTROL 
,• ISSUE LIMIT FORCE COMMAND 
,- ENABLE SERVOLOOP 


NOU UNDER ACTIVE CONTROL 


FLIH 

MBASE 

H,0 

MBASE+2 

A, 8 

MATH 

A, 6 

MATH 

A,TVO 

MBASE+4 

A,TV1 

MBASE+5 

A, TV2 

MBASE+4 

A, TVEXP 

MBASE+7 

A, 2 

MATH 

MBASE 

TRANS 

MBASE+2 

TRANS+2 


CALC TRANSITION VELOCITY 
» GET LIMIT FORCE CMD (FLI) 

;load into math board 

, -CONVERT TO FLOATING POINT 
,- SQUARE IT 

,’LOAD TRANSITION VELOCITY SCALE 
, -FACTOR INTO HATH BOARD 


,-MULT BY FLI**2 
J STORE TRANS VEL 
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3DEE 

3DEE 

3DEE 

3DEE 2A 88 3F 
3DF1 EB 
3DF2 3E 00 
3DF4 CD 84 3E 
3DF7 CD 99 3E 
3DFA 22 00 80 
3DFD 21 00 00 
3E00 22 02 80 
3E03 3E 08 
3E05 CD A5 3E 
3E08 2A 92 3F 
3E0B 22 04 80 
3E0E 2A 94 3F 
3E1 1 22 06 80 
3E1 4 3E 0A 
3E16 CD AS 3E 
3E1 9 DB A1 
3E1B E6 20 
3E1D CA EE 3D 
3E20 
3E20 
3E20 

3E20 2A 80 3F 
3E23 11 FD FF 
3E26 22 08 F7 
3E29 19 
3E2A DA 26 3E 
3E2D 21 00 00 
3E30 22 08 F7 
3E33 00 
3E34 C3 33 3E 
3E37 
3E37 
3E37 

3E37 21 00 00 
3E3A 22 08 F7 
3E3D 3E A1 
3E3F D3 EA 
3E41 3E 02 
3E43 CD 84 3E 
3E46 EB 
3E47 2A 8A 3F 
3E4A CD 99 3E 
3E4D D2 41 3E 
3E50 C3 10 3D 
3E53 
3E53 
3E53 


1870 * NOU HAVE TRANSITION VEL STORED AS FLOATING POINT, 32-BIT # 
1880 * SO START COMPARING THIS AGAINST (SINK RATE - U/G VEL) 

1890 * FOR DETERMINING START OF TRANSITION. 

; GET SINK RATE 
; S A VE IN DE 
; GET U/G VELOCITY 

; CALC : SINK RATE - U/G VEL 
, -CONVERT TO FLOATING POINT 


; ITS NOU IN MBASE+O, 1 ,2,3 
; LOAD TRANS VEL INTO HATH BOARD 


1900 L4 

LHLD 

SINK 

1910 

XCHG 


1920 

MVI 

A,MUX1 

1930 

CALL 

INI 

1940 

CALL 

SUB2 

1950 

SHLD 

MBASE 

1960 

LXI 

H , 0 

1970 

SHLD 

HBASE+2 

1980 

MVI 

A, 8 

1990 

CALL 

MATH 

2000 

LHLir 

TRANS 

2010 

SHLD 

MBASE+4 

2020 

LHLD 

TRANS+2 

2030 

SHLD 

HBASE+6 

2040 

MVI 

A,OAH 

2050 

CALL 

MATH 

2060 

IN 

MSTAT 

2070 

ANI 

20H 

2080 

JZ 

L4 

2090 * 



2100 * TRANSITION PHASE 

2110 * 



2120 

LHLD 

FLIM 

2130 

LXI 

D, -3 

2140 L5 

SHLD 

DACO 

2150 

DAD 

D 

2160 

JC 

L5 

2170 

LXI 

H,0 

2180 

SHLD 

DACO 

2190 L6 

NOP 


2200 

JMP 

L6 

2210 * 



2220 * TAKEOFF MODE 

2230 * 



2240 L2 

LXI 

H , 0 

2250 

SHLD 

DACO 

2260 

MVI 

A,0UT2 

2270 

OUT 

P0RT6 

2280 L7 

MVI 

A , MUX2 

2290 

CALL 

INI 

2300 

XCHG 


2310 

LHLD 

BXTHR 

2320 

CALL 

SUB2 

2330 

JNC 

L7 

2340 

JMP 

START 


; COMPARE AGAINST (SINK-U/G VEL) 

;GET STATUS, MASK 'LESS THAN' BIT 
; IS (SINK-U/G VEL) .LT. TRANS VEL? 
;N0, CONTINUE LOOPING 
YES, TIHE TO START TRANSITION 


;LOAD LIMIT FORCE COMMAND 

;SET RAMP RATE 

; OUTPUT CMD TO DAC 

; DECREASE LIMIT FORCE CMD 

;LOOP UNTIL CMD = 0 

;SET LIMIT FORCE CHD EXACTLY = 0 

; ST AY IN A LOOP UNTIL A RESET OCCURS 


; COMMAND A ZERO LIMIT FORCE 

JENABLE SERVO LOOP, LEAVE ENABLED 
; UNTIL STRUT POSITION LESS THAN THRESHOLD 
•GET STRUT POSITION 

;PUT IN DE 

;LOAD HL UITH THRESHOLD 
jCALCf STRUT - THRESHOLD 
: LOOP UNTIL STRUT EXTENDED FULLY 
;HAVE LIFTOFF, TURN OFF CONTROLLER 


2350 

2360 

2370 


* ROUTINE TO INPUT AND STORE DATA FROH 

* THREE MUX CHANNELS 
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3E53 




2380 * 




3E53 

3E 

01 


2390 IN3 

MV I 

A, MUXO 

J POINT MUX TO U/G ACCEL 

3E55 

21 

01 

F7 

2400 

LXI 

H,FCR 

POINT HL TO MUX/GAIN REGISTER 

3E58 

77 



2410 

MOV 

H,A 

JLOAD REGISTER 

3E59 

2B 



2420 

DCX 

H 

; POINT HL TO CHD/STATUS REGISTER 

3E5A 

36 

01 


2430 

M V I 

M, 01 

; START CONVERSION 

3E5C 

7E 



2440 HI 

MOV 

A, M 

,'READ STATUS 

3E5D 

07 



2450 

RLC 


,'DONET 

3E5E 

D2 

5C 

3E 

2460 

JNC 

Ml 

;NO f KEEP LOOPING 

3E61 

36 

00 


2470 

MVI 

M ,0 

; YES, RESET CONVERSION ENABLE 

3E63 

2A 

04 

F7 

2480 

LHLD 

ADDAT 

; GET DATA 

3 E 6 6 

22 

00 

80 

2490 

SHLH 

ACCEL 

; STO U/G ACCEL IN MATH BOARD, 

3E69 

22 

80 

3F 

2500 

SHLD 

FLIM 

; ALSO IN RAM 

3E6C 

3E 

00 


2510 

MVI 

A , MU X 1 

; REPEAT FOR U/G VELOCITY 

3E6E 

21 

01 

F7 

2520 

LXI 

H,FCR 


3E71 

77 



2530 

MOV 

M, A 


3E72 

2B 



2540 

DCX 

H 


3E73 

36 

01 


2550 

MVI 

H,01 


3E75 

7E 



2560 H2 

MOV 

A,M 


3E76 

07 



2570 

RLC 



3E77 

02 

75 

3E 

2580 

JNC 

M2 


3E7A 

36 

00 


2590 

MVI 

M 


3E7C 

2A 

04 

F7 

2600 

LHLD 

ADDAT 


3E7F 

22 

84 

3F 

2610 

SHLD 

UGVEL 

; STORE U/G VEL 

3E82 

3E 

02 


2620 

MVI 

A , MUX2 

; REPEAT FOR STRUT POSITION 

3E84 

21 

01 

F7 

2630 INI 

LXI 

H,FCR 


3E87 

77 



2640 

MOV 

M, A 


3E88 

2B 



2650 

DCX 

H 


3E89 

36 

01 


2660 

MVI 

H , 0 1 


3E8B 

7E 



2670 M3 

MOV 

A, M 


3E8C 

07 



2680 

RLC 



3E8D 

02 

8B 

3E 

2690 

JNC 

M3 


3E90 

36 

00 


2700 

MVI 

M,0 


3E92 

2A 

04 

F7 

2710 

LHLD 

ADDAT 


3E95 

22 

86 

3F 

2720 

SHLD 

STRUT 


3E98 

C9 



2730 

RET 



3E99 




2740 * 




3E99 




2750 * DOUBLE PRECISION 

SUBTRACT ROUTINE 

3E99 




2760 * HL 

= D E - H L 



3E99 

7B 



2770 SUB2 

MOV 

A,E 


3E9A 

95 



2780 

SUB 

L 


3E9B 

6F 



2790 

MOV 

L, A 


3E9C 

7A 



2800 

MOV 

A,D 


3E9B 

9C 



2810 

SBB 

H 


3E9E 

67 



2820 

MOV 

H,A 


3E9F 

C9 



2830 

RET 



3EA0 




2840 * 




3EA0 




2850 * ROUTINE TO SHIFT 

VALUE IN HL. LEFT 4 PLACES. 

3EA0 




2860 * 




3EA0 

29 



2870 SHL 

DAD 

H 


3EA1 

29 



2880 

DAD 

H 
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3EA2 

29 



2890 


BAD 

H 



3EA3 

29 



2900 


BAD 

H 



3EA4 

C9 



2910 


RET 




3EA5 




2920 

♦ 




. 

3EA5 




2930 

* ROUTINE TC 

1 ACTIVATE HATH BOARD AND WAIT FOR RESULT 

3EA5 




2940 

* ACCUH HAS 

OPCODE. 



3EAS 




2950 

* 





3EA5 

D3 

AO 


2960 

HATH 

OUT 

IOBAS 

; COMMAND HATH BOARD 

TO START 

3EA7 

DB 

A7 


2970 

UA1T 

IN 

IOBAS+7 

; GET FLAG BYTE 


3EA9 

E 6 

01 


2980 


ANI 

01 

;CHECK BUSY BIT 


3EA6 

C2 

A7 

3E 

2990 


JNZ 

WAIT 

;STAY IN LOOP UNTIL 

NOT BUSY 

3EAE 

C9 



3000 


RET 




3EAF 




3010 

* 





3EAF 




3020 

* 





3EAF 




3030 

* SPECIAL CHECK-OUT 

ROUTINES 


3EAF 




3040 

* 





3EAF 




3050 

* 

ROUTINE 

TO INPUT 

' A VALUE FROH A/D, 

STORE IN RAH. 

3EAF 




3060 

* 





3EAF 

F3 



3070 


BI 




3EB0 

3E 

00 


3080 


MV I 

A, 00 

; SELECT CHAN 0 


3EB2 

CD 

84 

3E 

3090 


CALL 

INI 



3EB5 

CF 



3100 


RST 

1 



3EB6 

00 



3110 


NOP 




3EB7 

00 



3120 


NOP 




3EB8 




3130 

* ROUTINE TO DO PGA 

TEST ON A/D 


3EB8 

F3 



3140 


BI 




3EB9 

21 

01 

F7 

3150 


LXI 

H,FCR 



3EBC 

36 

00 


3160 

PGA 

M VI 

H , OO 



3EBE 

36 

CO 


3170 


HVI 

M r OCOH 



3EC0 

C3 

BC 

3E 

3180 


JNP 

PGA 



3EC3 

00 



3190 


NOP 




3EC4 




3200 

* 

ROUTINE 

TO OUTPUT A VALUE TO DACO, 

BACK 

3EC4 

F3 



3210 


BI 




3EC5 

00 



3220 


NOP 




3EC6 

21 

00 

00 

3230 

R2 

LXI 

H,0 



3EC9 

22 

08 

F7 

3240 


SHLD 

DACO 



3ECC 

22 

OA 

F7 

3250 


SHLD 

DAC1 



3ECF 

00 



3260 


NOP 




3ED0 

00 



3270 


NOP 




3ED1 

00 



3280 


NOP 




3ED2 

C3 

C6 

3E 

3290 


JHP 

R2 




? 
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Appendix B 

ELECTRONIC CONTROLLER DETAILED DESCRIPTION 
GENERAL CONSTRUCTION 

The controller consists of four boards, together with + and - 5V dc power 
supplies, control switches, and test jacks. The four boards perform the 
following functions. 

(1) Central Processing Unit (CPU) Board: The CPU board is part 
of the computer and performs the basic digital computations as well as the 
logical computations which determine the operating mode of the controller. 

(2) Arithmetic Board: The arithmetic board performs, in digital 
form, the required multiplication and division functions associated with 
energy and transition velocity. 

(3) Linear (Analog) Board : The linear board provides the control 
laws and functions associated with the force and position loops. It also 
computes the wing/ gear velocity by integrating the wing/ gear acceleration. 

In addition, it incorporates the switching circuitry which is actuated by 
signals from the CPU. 

(4) Analog-to-Digital (A/D) and Digital -to -Analog (D/A) Board ; 

This controller board converts the sensed and -computed analog quantities 

to digital form so that the digital computations can be performed to determine 
the commanded limit force, and converts the limit force to analog form so 
that it can be used as the input to the force loop. 

The physical location of the units is shown in Figure B-l and the 
functional interrelationship (block diagram form) is shown in Figure B-2. 

As received from the factory, the A/D board is configured to operate 
over the range of 10 to +10 V. All of the analog signals in the controller lie 
between 0 and +10 V. Therefore, to achieve maximum accuracy, the A/D 
board was reconfigured to operate in this range in accordance with the manual's 
jumpering instructions. The board is sent to NASA in this configuration. 
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Analog Board 



Figure B-l. - Landing Gear Controller (Rear View) 





cc 
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Figure B-2, - Controller Static Gains and Scale Factors 







CONTROL LAWS 


The control laws implemented in the controller are presented below in 
Figure B- 3. {Refer to SYMBOLS, Page 3, Report Proper.) 



Figure B-3. Controller Control Laws 


As shown in this figure, the force loop compensation consists of a notch 
filter with a center frequency of 251. 3 rad/sec and two lead-lag networks, 
0.0281S + 1/0.0141S + 1 and 0.001S + 1/0.0001S + 1. The position loop, 
which insures that the strut returns to its static position, incorporates a 
simple lag network: Kf/O.lS+1. 

GAINS AND SC A LING -LIN EAR 

The static gains and scale factors of the controller are shown in Figure 
B-2. These gains and scale factors were chosen to be compatible with the 
maximum values of the system parameters and the 10~V maximum of the 
microprocessor. 
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GAINS AND SCALING - DIGITAL 


The digital scaling is accomplished as follows: 

(1) W/G acceleration (X wg ) 

x wg " 6 (1. 8) V/g • 409. 5 bit/V = 4423 bit/g, or 0. 000226 g/bit 

(2) W/G velocity (X wg ) 

X^ g = 3. 937 V/m/ sec • 409. 5 bit/V = 1612 bit/m/ sec (40. 95 bit/ in/ sec) 
or 6.203 (10~ 4 ) m/ sec/bit (0.02442 in/ sec/bit) 

(3) Strut displacement (X s ) 

X s = 0. 715 * 39. 37 V/m* 409. 5 bit/V = 1. 153 (10 4 ) bit/m (292.8 bit/ in) 
or 8. 673 (10” 5 ) m/bit (0. 003415 in/bit) 

(4) Strut velocity (V s ) 

V s = 3.937 V/m/ sec* 409. 5 bit/V = 1612 bit/m/ sec (40. 95 bit/ in/ sec) 
or 6.203 (10" 4 ) m/ sec/bit (0.02442 in/ sec/bit) 


(5) 


Work potential of the strut (WP) 

WP - F wg • <X smax -X s ) . Mx wg (X smax -X s ) 

If X W g = lg and X s - 0. 0254 m (1 in) 

WP = Mg = W = 363. 3 N • m (3215 in. lb) 

In digital terms, 

WP = 6 (1.8) V* 409. 5 bit/V* 0. 715V • 409. 5 bit/V 
= 1.2949 (10‘ 6 ) bits 


Therefore, the scale factor of WP is: 

1.2949 (10 6 ) bits 

^63. 3 N • m 3564 blt / N ’ m; or 2.806 (10 -4 ) N • m/bit 


(402. 7 bit/ in • lb, or 0. 002483 in • lb/bit) 
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(6) Kinetic energy (KE) _t 

KE = 1/2 • W/ g (VrpQrp)^ where VrpQip - X W g touchdown \ ^ wg 
If V = 0.0254 m/ sec (1 in/ sec) ° 

KE = 0. 4706 N • m (4. 1645 in* lb) 

In digital terms, 

KE = 0. 1 V (409. 5 bit /V) 2 = 1676. 9 bits 

Therefore, the scale factor of KE is 

1676. 9/0. 4706 = 3564 bit/N • m, or 2. 806 (10 -4 ) N • m/bit 
(402. 7 bit/ in • lb, or 0.002483 in • lb/bit) 

Which is the same scale factor as that for WP, and the 
two terms can be compared directly. 

(7) Decrease of limit force command (Fyi) during the transition 
from impact phase to rollout phase: 

The scale factor of F^j is 1. 324 (10 4 )N (2977 lb) for 10 V; and, 10 V 
corresponds to 4095 digital bits. Therefore, the digital 
scale factor for Fyj is 4095/1. 324 (10 4 )N; or 0. 3094 bit/N 
(1. 376 bit/ lb). During transition, Fy is decreased at 

a rate of 1. 379 (10^) N/sec (31 000 lb/ sec), or digitally at 
1. 379 (10 5 ) 0. 3094 = 42 642 bit/ sec. 

(8) Transition Velocity: 

From Figure 3.2. 3 of the system specification: 

v - F LI 2 
T " 2 (W/g) R 

where W = aircraft weight per gear and R is the limit 
force transition rate. 


The scale factor of V T is determined as follows: 

2 , 

W/g = 1459 N * sec /m 

Then 1 Newton of F T T produces 
2 ^ 

(1) - o a oa tv* / o or* t 


2 (1459) (1.379 * 10 ) 


= 2.486 (10 ) m/sec of V r 
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Digitally, the scale factor for F LI (from the previous section) is 0, 3094 
bit/N (1. 3755 bit/lb) Then, the scale factor for V t /F li = 

(0. 3094) 2 = 0. 0956 bit/N (1. 892 bit/lb) of F LI 
Therefore, the scale factor for V-p is: 

2"q = 3. 845 (10 7 ) bit/m/sec (977 000 bit/in/sec) 

2. 486(10 y ) 


This scaling must be matched to the scaling of V TOT ; that is. 


/ 


X wg dt + V g . 


0. IV 409. 5 bit 

The scalmg of Vtat = n “X oTZ 7 x 77 

6 1UT 0.0254 m/sec V 

- 1612 bit/m/sec (40. 95 bit/in/sec) 

To provide this scaling for V-p it must be multiplied by 1612/3. 845 (10 7 ) 
= 0.00004191 using the arithmetic board. This is accomplished as follows: 


0.00004191 DECIMAL (D) = 

.00000000000000101011111100100010000101 BINARY (B) 
or 1. 01011111100100010000101 x 2 -15 <B) 

The exponent is -15 (D) 

The bias in the arithmetic board is 

07F HEXIDECIMAL (H) or 127 (D) 

Therefore, the number must be applied to the arithmetic board with a 
bias of 127-15 = 112 (D). In addition, a factor must be applied due to the 
fact that the numbers from the A/D converter are stored in the most signi- 
ficant 12 bits out of 16 so that the number for V TO t is a factor of 16 too high. 
The transition velocity (V-p) is a function of Flj 2 # and is a factor of (16) 2 too 
high. The net result is that the number for V T is too high by a factor of 16. 

It must therefore be reduced by a factor of 16 or 2^. 

Therefore, the exponent of the applied number is 112-4 or 108 (D) = 6C (H) 
= 01101100 (B). 
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A sign bit ("0" for positive) must precede the exponent. The format 


of the applied number, is: 

Sign Exp. Fraction 

oouon^ 

M+3 M+2 M+l M 

Therefore, if this number starts at location M, 


the contents of memory 


are: 


M 85(H) 

M+l C8(H) 

M+2 2F(H) 

M+3 36(H) 

LINEAR CIRCUIT DESCRIPTION 


The linear circuit is shown in HR drawing 88000080-201. Power for the 
linear components is obtained from two auxiliary power supplies which 
provide +15 Vdc and -15 Vdc. The signal from the wing / gear accelerometer 
is applied to a differential pair of amplifiers, U22 and U23. This approach 
was taken in order to minimize the noise on the low-level signal. The output 
of the differential pair is applied to U4 and the output of U4 is biased by 
potentiometer R6 to provide a 1-g offset level. The biased signal is applied 
to U3 and then to U5 which has a gain of 6 and the output of which is the 
wing/ gear interface force which is applied to the microprocessor. The 
acceleration signal is available at J13 on the front panel. 

The acceleration signal from U3 is also applied through potentiometer 
R2 to integrator U2, the output of which is the wing/ gear velocity signal and 
is applied to the microprocessor. R2 provides a means of adjusting the inte- 
grator gain. The integrator is enabled by analog switch U1A which removes 
the short circuit across the capacitor upon receiving an enable signal from 
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the microprocessor. The wing/ gear velocity signal is available at J12 on the 
front panel. Switch S7 is provided on the front panel in order to allow a 
simulated wing/ gear signal to be applied to the microprocessor for test pur- 
poses. The simulated test signal is applied to J25 on the front panel. 

The commanded limit force is algebraically summed with the wing/ gear 
force by means of R22 and R2 3 to produce the force error which is then 
amplified by U6 and U7. The limit force command signal is available at J15 
on the front panel. Switch 5 allows a simulated limit force command signal 
(applied at J23 on the front panel) to be used for test purposes. R15, in the 
feedback path of U7, allows the forward loop gain of the system to be adjusted 
as required. 

The output of U7 is applied to the notch network (bridged T) which is 
composed of R16, R17, C2 and C3, and the output of which is applied to U8. 

The output of U8 is applied to U9, which provides one of the lead-lag functions, 
and then through U10 to Ull, which provides the other lead-lag function, and 
the output of which is applied to U12. The signal from U12 passes through 
one path of dual analog switch U1 and then to U14, the output of which is the 
servovalve command signal, and is applied to the servovalve controller. This 
signal is available at J19. Ul closes the forward loop path upon receipt of an 
enable signal from the microprocessor. R31 on the front panel provides a 
means of biasing the servovalve. 

The strut position is set by R41 on the front panel. The signal from this 
potentiometer passes through U17 and U16 (when the servoloop is enabled) 
and is algebraically summed, at U15, with the signal from the strut poten- 
tiometer, after it has passed through U19 and U18. The output of U15 is the 
strut position loop error. It passes through U21 and is applied to the force 
loop at U6. The strut position command signal is available at TP1, the strut 
position signal is available at J18 and the strut position error signal is 
available at J17, all on the front panel. The controller is enabled and reset 
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by means of switches on the front panel which are provided for this purpose. 

In addition, the strut hydraulic pressure signal and pneumatic pressure signal 
are available at J21 and J22 respectively on the front panel. 

The signals for the controller are applied at the rear of the unit. These 

are: 

J2 - 28 Vdc 

J4 - W/G acceleration 

J5 - Hub acceleration 

J6 - Strut position 

J7 - Strut pneumatic pressure 

J8 - Strut hydraulic pressure 

J9 - Servovalve command 

In order to set and maintain the initial hydraulic pressure in the gear an 
auxiliary pressure loop is used prior to enabling of the servoloop. To 
accomplish this, the pressure signal is amplified by amplifiers U24, U2 5, 
U28, U30 and TJ32, the output of which is added to the servoloop command 
signal through switch U33. When the servoloop is disabled, the switch is 
closed and allows the pressure signal to close the loop. The pressure is then 
controlled by a servovalve bias signal. When the servoloop is enabled the 
switch is opened and the pressure is free to vary in response to the loop 
command signal. 


DIGITAL SOFTWARE 


The digital software program is listed in Appendix B. Appended to this 
listing are routines for testing the arithmetic board and A/D board. 

The program is in the C. P. U. twice; that is, the C. P. U. contains two 
PROM's, each containing the entire program. One PROM is at location 0000 
and is the one normally used. No special procedures are required to use it. 
When power is applied the computer starts at this location, and once the 
Controller Enable signal is received it assumes control of the process. 
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The second PROM is intended for test and program changes if required. 
It is located at address 0800 in program memory but is programmed to start 
at address 3D 10 in RAM. To use it, the first PROM must be replaced by the 
monitor ROM (at 0000) and then the program can be controlled by a standard 
teletypewriter connected to the proper socket on the rear of the CPU board. 
The contents of locations 0800 to 09F9 are moved to new locations starting at 
3D 10 with a teletype input: M0800, 09F9, 3D10 RETURN. Then, any input 
desired can be applied to the computer by means of the teletype -- for testing 
or for program changes. To operate in this mode an input is required -- 
G3D10 RETURN -- before any test. 

If permanent program changes are required the PROM must be "burned" 
to contain the new program. . 
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Appendix C 
TEST PROCEDURE 

GEAR CHARGING PROCEDURE 

(1) With the gear vertical, and the dead weight of the beam as a static 
load, bleed any accumulated gas from the hydraulic port of the gear until 

hydraulic fluid escapes from the port. 

(2) Bleed gas or hydraulic fluid from the pneumatic charge port of 

the gear until the gear is fully compressed. 

(3) Recheck the hydraulic port for any additional accumulated gas. 

(4) If hydraulic fluid does not emerge from both gear ports in the 
fully compressed condition then fluid must be added. One method of 

accomplishing this is as follows: 

(a) Turn on the controller 28-Vdc supply and electronics. 

(b) Press the controller "RESET" button. This applies 
a positive bias command (pressure bias) to the servovalve. 

(c) Turn on the gear hydraulic supply pump. Momentarily 
raise the pressure by means of the main relief valve to approximate- 
ly 4140 kPa (600 lb/in^), and then reduce this pressure to about 

690 kPa (100 lb/ in 2 ). 

(d) Slowly open the gear isolation valve (see Figure C-l). 

This should apply hydraulic pressure to the gear and bleed gas at 
both the hydraulic and pneumatic ports of the gear until hydraulic 

fluid escapes from both. 

(e) Close the gear isolation valve. 

(f) Turn off the gear hydraulic supply pump. 

(g) If the gear has extended during this procedure repeat 

steps 1, 2 and 3. 
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(5) Connect a nitrogen charging system to the pneumatic charge 
port and apply the desired pre-charge pressure. 

NOTE: The gear has a working pressure rating of 3450 kPa (500 lb/in ). 

Therefore, caution should be exercised to avoid exceeding the 

proof pressure of 5170 kPa (750 lb/ in ). 

(6) Until the strut is fully extended and is stabilized at the 
desired charge pressure, slowly lift the upper gear using the load system. 

(7) Close the pneumatic charge port and remove the charging 
system. 


PASSIVE GEAR TEST PROCEDURE 

(1) Turn on the load system electronics and allow 30 minutes 
for warm-up. 

(2) Turn on the load hydraulic supply system. Be sure minimum 
capacities are 114 liter/min (30 gal/min) and 2. 413(10 4 ) kPa (3500 lb/ in ). 

(3) Move the mode switch on the load controller to "POSITION . 

(4) Open the load system supply valve and slowly raise the load 
system pressure to 1. 45 (10 4 ) kPa (2100 lb/in 2 ). If the position command 
potentiometer has been preset to maximum drop height, the beam will move 
to its retract stop, thereby raising the gear. 

(5) If the position command potentiometer has not been pre-set, 
then slowly adjust it to position the beam against the retract stop and 
check to see if the gear charge pressure is at the desired value. 

(6) Set the recorder to monitor the required parameters and set 

the channel gains. 

(7) Momentarily move the reset/ operate switch on the load 

controller to "RESET" and then return it to "OPERATE". 

d 2 

(8) Raise the load system pressure to 2. 069 (10 4 ) kPa (3000 lb/ in ). 
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(9) Start the recorder and move the mode switch to "VELOCITY/ 
LOAD" to drop the gear. 

(10) After the drop, reduce the load system pressure to 2100 
lb/ in 2 and move the mode switch to "POSITION". 

(11) For further testing repeat steps 6 through 9. 

(12) When testing is concluded, reduce the load system pressure 
to minimum and close the load system supply valve. 

ACTIVE GEAR TEST PROCEDURE 


(1) Turn on all load system, controller, and servovalve 
controller electronics. Allow 30 minutes for warm-up. 

(2) Follow the procedures of steps 2 through 7 of "Passive 
Gear Test Procedure". 

(3) Place the gear controller in the "RESET" state. 

(4) With the gear isolation valve closed, turn on the gear 
hydraulic supply pump and adjust its main relief valve to provide 
6895 kPa (1000 lb/ in ) as read on the supply pressure gage. 

(5) Adjust the "BIAS" control on the servovalve controller to 
produce approximately the same gear supply pressure as the gear 
charge pressure. 

(6) Slowly open the gear isolation valve while observing the 
pressure in the gear. The gear internal pressure must not fall more 

than 345 kPa (50 lb/ in 2 ) below its pre-set charge pressure. Otherwise, the 
gas may be forced into the hydraulic side of the gear and necessitate re- 
charging m accordance with the first section, "Gear Charging Procedure", 

(7) Readjust the "BIAS" control on the servovalve controller to 
obtain the desired gear internal charge pressure, as read by the gear 
hydraulic pressure transducer. 
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(8) Confirm that the "S. V. CMD. BIAS" control on the gear 
controller is set to the proper value for the gear charge pressure. 

(9) Set the sink speed value by means of the "SINK SPEED control. 

(10) Set the load "VELOCITY" command potentiometer for the 
desired sink speed. 

(11) Momentarily move the reset operate switch on the load 
controller to "RESET" and then return it to "OPERATE". 

(12) Raise the gear supply pressure to 2.069 (10 4 ) kPa (3000 lb/in ) 
and recheck the gear internal pressure. If necessary, readjust the servo- 
valve controller "BIAS" control for the desired gear pressure. 

NOTE: To prevent overheating do not operate the gear 
hydraulic supply at 2. 069 (10 4 ) kPa (3000 lb/ in ) until 
immediately prior to the drop. 

2 

(13) Raise the load system pressure to 2.069 kPa (3000 lb/ in ). 

(14) Start the recorder and move the load mode switch to 
"VELOCITY LOAD" to drop the gear. 

(15) After the drop, reduce the gear supply pressure to 

6900 kPa (1000 lb/ in 2 ). 

(16) Reduce the load system pressure to 1. 45 (10 4 ) kPa (2100 lb/ in ). 

(17) Close the gear isolation valve. 

(18) Raise the gear by moving the mode switch on the load 
controller to "POSITION". 

(19) Return the gear controller to the "RESET" mode. 

(20) Slowly open the gear isolation valve. The gear internal 
pressure should return to the pre-set charge value. For further 
testing repeat steps 9 through 19. Otherwise proceed to step 21. 

(21) Reduce the load pressure slowly to minimum allowing 

the gear to settle gently, and close the load system supply valve. 
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(22) Turn off all electronics except the recorder. Wait about 
two minutes (to allow capacitors to discharge) and run a short record 
for channel zero references. 



Appendix D 

OSCILLOGRAPH RECORDINGS 
Ke y to Recordings 


Channel 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Parameter 


Sensitivity 


Zero Position 
(with respect 
to reference line) 


W/G Velocity 

Net Force 
(Accelerometer #1) 

Servovalve Spool Pos* 

Limit Force Command 

Net Force (filtereo 
(Accelerometer #2) 

Strut Position 
Servovalve Command 
Gear Hydraulic Pressure 
Gear Pneumatic Pressure 
Lift Force Command 
Lift Simulator Position 
Lift Force 


0*25 m/sec/cm (0.25 in/ sec/ in) 
2224 N/cm (1270 lb/ in) 

0 . 05 cm/ cm (0. 05 in/ in) 

876 N/cm (500 lb/ in) 

2224 N/cm (1270 lb/ in) 

0. 0127 m/in (1.27 in/ in) 

1, 97 ma/cm (5 ma/in) 

277 kPa/cm (102 lb/in 2 in) 

2 77 kPa/cm (102 lb/ in 2 /in) 

1. 97 V/ cm (5 V/in) 

0. 05 ml cm (5 in/ in) 

8756 N/cm (5000 lb/ in) 


16. 76 cm (6. 6 in) 

7. 62 cm (3. 0 in) 

15. 24 cm (6. 0 in) 

2. 54 cm (1. 01 in) 
10. 19 cm (4. 01 in) 

19. 71 cm (7. 76 in) 
9. 17 cm (3. 61 in) 
19. 81 cm (7. 8 in) 
19. 30 cm (7. 6 in) 
18. 26 cm (7. 19 in) 
12. 70 cm (5. 0 in) 
17. 78 cm (7. 0 in) 


148 



Reference Line 
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. 1048 kPa (152 lb/ In 2 ), 1. 524 m/sec (5 ft/sec). Active 
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